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development and application of a non-gaussian atmospheric 
turbulence model for use in flight simulators 
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The simulator experiment produced results which were not 
conclusive. Although the pilots judged that the simulator 
response to the turbulence generated by this new method indi- 
cated that this model produced turbulence somewhat more real- 
istic ^ patchy", and less monotonous than that generated from 
the Gaussian model, the other subjective information obtained 
as pilot opinion was inconclusive regarding distinguishing 
the various types of turbulence inputs. It is concluded that 
more extensive experience is required to evaluate the worth 
of this new model in handling qualities studies. 


INTRODUCTION 


The widespread use of flight simulators to evaluate 
vehicle handling qualities, pilot workload, ride quality, 
and other design factors which must be evaluated with a human 
in the control loop, has resulted in the development of 
highly sophisticated simulators, including six degree-of- 
freedom moving base simulators and complex variable stability 
aircraft. ^ 


It has been found that realism requires the presence of 
an external disturbance representing atmospheric turbulence. 

A ^®^listic representation of turbulence becomes especially 
important in the simulation of STOL vehicles because of their 
sensitivity to turbulence, as even light to moderate turbu- 
lence may seriously degrade the controllability and ride 
quality of these vehicles. 

_ methods of generating turbulence signals for use 

in simulators have been developed. Some of these will be men- 
below, along with their shortcomings. 

The most direct method is the use of recorded gust time 
histories obtained from flight experiments (ref. 1) . This 
form is certainly realistic, and may even be especially use- 
ful when the turbulence is generated by terrain and obstacles, 
but for general studies it has several disadvantages. It is 
difficult to modify such records to fit anything except the 
specific situation recorded, which may or may not be typical, 
there is a possibility that pilots may become familiar 
wfth the time histories if they are used repeatedly. There- 
fore a large number of recordings representing a wide variety 
of turbulence conditions becomes necessary. 

A more useful turbulence simulation might result if one 
could generate non-repeating random time histories with pro- 
per statistical properties of the actual atmospheric turbu- 
lence being simulated. Such a method should allow continuous 
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adjustment of these properties in order to simulate changes in 
the flight environment. Three such techniques are described. 

The first of these represents turbulence by a summation 
of sinusoidal waves with various amplitudes and frequencies 
(ref. 2) . The principal objection to this method is that it 
attempts to model a process containing an infinite nvimber of 
frequencies by using a small number of sinusoids, and there- 
fore valuable information may be lost. Nevertheless, the sum 
of sine waves model does produce a non-repeating time history 
and has been used in flight simulation. 

A second method generating random time histories uses a 
linearly filtered Gaussian white noise signal (ref. 3). This 
Gaussian model or Gaussian simulation requires only inexpen- 
sive equipment and is not particularly complicated if the 
gust components are uncorrelated. This Gaussian model has 
two major deficiencies. First, the artificial turbulence 
produced is necessarily Gaussian, whereas real turbulence 
is non— Gaussian . This results in an underestimation of the 
large gusts by the Gaussian model which is a serious defi- 
ciency in that the large gusts cause the greatest control 
problems and result in the greatest passenger discomfort. 
Secondly, the Gaussian model does not truly represent the 
patchy nature of real turbulence. Turbulence occurs in 
patches of relatively intense activity separated by intervals 
of relative calm. These deficiencies result in a model which 
produces turbulence inputs which pilots classify as too con- 
tinuous or monotonous compared to real turbulence. 

A third method, proposed by one of the authors of this 
paper (ref. 4) consists of using the product of two Gaussian 
processes, one representing a constant intensity turbulence 
time history and the other representing the time varying stan- 
dard deviation. This method produces a non-Gaussian model and 
has the same versatility as the Gaussian model. However, the 
simulation is not as realistic as might be hoped. The approach 
described in this paper is a modification and extension of 
this third method. 

The purpose of the study described in this report was to 
develop a technique for producing artificial turbulence which 
would overcome the deficiencies of the techniques described, 
and yet retain the practical advantages of the Gaussian simu- 
lation in ease of application. Such a technique was devel- 
oped and put into a form which could be used in simulators. 

The application of this new non-Gaussian model is somewhat 
more complex than that of the Gaussian model, but the new 
model does include the possibility of producing the large 
gusts and patchy characteristics of turbulence lacking in the 
Gaussian model. 
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This new model and others were mechanized to be included 
in a flight simulator experiment in order to determine to what 
extent pilots are sensitive to changes in the "patchiness" of 
a turbulence model. Pilots were asked to fly a constant alti- 
tude tracking task for a typical STOL aircraft in the landinq 
configuration with four different turbulence models, recorded 
real atmospheric turbulence, a Gaussian model, and two non- 
Gaussian models of differing patchiness. These models differed 
only xn thexr probability densities. 

The following sections of this report describe the pro- 
posed turbulence simulation and the flight simulator experiment 
xn detaxl. Brxefly the sections consist of: 

a.) Description of Turbulence Simulation. This section 
describes the characteristics of real turbulence and several 
sxmulatxons now in use. The proposed new simulation is then 
dxscussed and it is shown that this simulation produces an 
improved representation of atmospheric turbulence. The manner 
in which the statistics of the spatial distribution of gusts 
over the vehicle surface are obtained is presented. The 
section also includes the method used to introduce the artifi- 
cial turbulence into the flight simulator. 

4. Simulator Experiment. Here are found details 

ot the flight simulator experiment, in which pilots were asked 
to rate the realism of the various turbulence models. 

c. ) Results. This section summarizes and discusses the 
results of the previous sections. 

d. ) Conclusions. 


e.) Appendixes. Five appendices are included. 

Appendix A describes the details of the turbulence 
simulations used in the flight simulator experiment. 

Appendix B describes the equations of motion used 
in the experiment. 

Appendix C describes the flight simulator motion 
compensation. 


Appendix D presents the pilot question sheet. 

Appendix E presents additional results of the pilot 
questionnaire not presented in the main text. 
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SYMBOLS 


a 

b 

c 

C, 


CT 


CT 


random function of time 

wingspan, m ; also, a random function of time 

wing chord, m ; also, a random function of time 

= Drag — ^ nondimensional drag coefficient of the 

(1/2) pu^^S 
airplane 

, 3C 

_ Lift 

(1/2) pu^^S 
airplane 

= f— ") 


nondimensional lift coefficient of the 


T' 0 


= Moment ^ nondimensional rolling moment 


(1/2) pu^^Sb 


coefficient of the airplane 


,2u„ ,9C„ 

( o\( £ 


V b 


9r ) o 



contribution 

to 

from 

the 

distribution 

of head-on 

ug 

gusts across 

the 

wingspan 





contribution 

to 

from 

the 

distribution 

of vertical 

wg 

gusts across 

the 

wingspan 
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(^) 

r 0 

Pitching Moment 

5 , nondimensional pxtching-moment 

(1/2) PUq^Sc 

coefficient of the airplane 


^^m^ 



1 

0 


9C , 

> c / 

\ 9q 

2 



( \ — )• (pitching moment about airplane c.g. due 

PUq Sc 

to horizontal tail only) 

/ 9C . . 

<"o> ^-5?)o 


(^)(^) 

9a 0 

/ 9C V 

I— ^ 

\ 9a Jo 

9a 

/ 9C V 

(— ^ 

'96e / 0 

Yawing Moment 


(l/2)puQ^Sb 

coefficient of the airplane 
f 2u.\ / 9C X 
' b A 9p Jo 


nondimensional yawing-moment 
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2u_. /3C 


/ \ / 


pu ^Sb 

0 


. (Yawing moment about airplane c.g. due 


to vertical tail only 


/3C . 

' nt 


V 83 


contribution to C from the distribution of head-on 

n 

gusts across the wingspan 

contribution to C from the distribution of vertical 

n 

gusts across the wingspan 
/8C V 
^ 83 ^ 


,'8C 

I— 

^ 8r 


Thrust 
(1/2) pu^^S 

^8 6t /o 


, thrust coefficient 


autocorrelation of random process u(t) 

cross correlation of random processes u(t) and v(t) 


PUq^S 


. (x-component of aerodynamic force acting 


on the airplane) 
8C , 

i — ^ 

\ 8a jo 


y — • (y-component of aerodynamic force acting 
^Pu 

on £he airplane) 
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2u^ , 3C 

ON'' V \ 




b A 9p A 


AC 


V b )\ dr Jo 


t3 


^pu 

0 

/9C ^ 

_ '3C ^ 

^T?)o 


(y-force due to vertical tail only) 


6a 


6r 


,3C 

^8 6a ^0 


8C 


( u 

96r 


^2 ( T~) * (z-component of aerodynamic force acting 


2^ 
PUq S 


on the airplane) 


= ( 


2u^. , 3C . 
o\ / z\ 

c J\ 3q Jo 


= i 2 N 


'ta 


'ta 


a 


a 


[ 2“j * due to horizontal tail only) 


pu S 

0 

/ 9C ^ 

( zt \ 

V 9a jo 


V c J\ • j ^ 
' 9a 0 

, 9C 
{ z^ 

\ 9a j 0 
2u„, , 9C 




Z\ 


9a ^ 0 
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.9C , 

= ( 

Vs 6 e )o 


random function of time 
frequency, hertz 

p ^ p ) X, y, and z-components of engine force, N 

' E E 

y z 

p ^ p ) X, y, and z-components of total force on 

"^y ' '^z the airplane (excluding gravity) , N 

maximum engine power, W 

modified Bessel function of second type, order zero 

longitudinal scale-length of atmospheric turbulence, m 

distance from airplane c.g. to tail mean-aerodynamic- 
center, m 

airplane mass, kg 

^ ]yi ) X, y, and z-components of moment due to 
^z engines, kgm 

normalized central moment 

^24 ) x,y, and z-components of total moment on 

"^y '^z the airplane (excluding gravity) kgm 

integer 

roll rate (right wing down, positive) , rad/sec 
probability density 

cumulative probability distribution 
pitch rate (nose-up, positive) , rad/sec 
= (1/2) , dynamic pressure, N/m^ 

yaw rate (nose right, positive) , rad/sec 
Laplace transform variable 
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C • 2 

o wing area, m 

Sears function for a two-dimensional wing in a 
sinusoidal upwash 

STOL short take-off and landing 

t time, sec 


T total time history length (sec) 

Uq reference flight speed in still air, m /sec 

(u, V, w) X, y, and z-components of airplane velocity with 
respect to still air, m /sec 


<“b' ''b' 


"b> 


Y' 9-nd X— components of airplane velocity with 
respect to gusty air (i.e., u_ = u + u , etc ) 
m /sec ° 


(u 


g 


C 



w 


X, y, and z-components of gusts at the 
^i^pl^ne c.g. (u , v , w are positive 

^c ^c 

when in the direction of negative x, y, and 
z) , m /sec 


^^g/ ) components of gusts at the airplane tail 

^ ^ t (same sign convention as c.g. gusts) , m /sec 

^e' north, east, and down components of the airplane 

velocity with respect to still air, m /sec 


/ 2 2 2 

^RW ■*" ''b ' airplane speed with respect to 

gusty air, m /sec 

(x, y, z) distances measured in airplane body axes; (body— axes 

are centered at aircraft c.g., attached to the airplane, 
and have x-axis through the nose, y-axis out the right 
side of the airplane, and z-axis downward with respect 
to the airplane) , m 

distance airplane has travelled in the north, 
east, and downward directions, m 


a 


_1 

- tan — , airplane angle of attack, rad 
B 
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B 

6 

a 

6 

e 

6 

r 

6T 

n 


angle o£ attack o£ the zero-lift line of the wing, rad 
V 


tan 


-1 


B 




sideslip angle, rad 


angle moved by one aileron surface, rad 
(positive 6^ gives positive rolling moment) 

angle moved by elevator surface (trailing-edge up, 
positive) , rad 

angle moved by rudder surface (trailing-edge right, 
positive) , rad 

pilot's throttle movement (throttle forward is 
positive) , measured in % of full throttle 

Gaussian white noise 


a 

P 




standard deviation 

3 

atmospheric density, kg/m 


power spectral density of the variable i , m 
(units of the variable i) /hertz 


in 


cross spectral density ^hertz 

(units of variable i) (units of variable 3 )/hert 


Subscripts: 

u, V, w denotes u, v, or w gust component 


= value at reference-flight condition with no wind 
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TURBULENCE SIMULATIONS 


In this section the general problem of statistically 
modeling atmospheric turbulence is discussed, and detailed 
descriptions are given of various techniques for accomplish- 
ing this modeling, including the proposed non-Gaussian model. 

First, basic definitions of statistical properties of 
atmospheric turbulence are reviewed. Then the statistical 
properties of atmospheric turbulence which are to be used in 
the various models are discussed in the light of presently 
available information. Next, presently used simulation tech- 
niques are discussed in some detail. Finally, a non-Gaussian 
simulation method is proposed, and it is shown that this tech- 
nique can closely model the statistical properties of real 

including spatial effects, providing proper informa- 
tion and data regarding the real turbulence are available. 


Review of Basic Definitions 

The following is a collection of definitions applying 
to the statistical quantities used throughout this report. 

Stationarity . - A random process is stationary if its 
statistical properties are not dependent on the time of their 
measurement. One could, for example, collect an infinite 
number of time histories which are representative of the pro- 
cess. This set of histories is called an ensemble, and could 
look as below. 
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If one takes an ensemble average across the ensemble, 
at tj^ and at ^2 / etc., and if these ensemble averages 

are not functions of time the process is stationary. 

All of the definitions which follow assume stationarity , 
with the exception of homogeneity. 

Homogeneity . - A random process is homogeneous if its 
statistical properties are independent of position. Thus, if 
atmospheric turbulence were homogeneous (which it is not) , 
its statistical properties would be independent of altitude 
and geographical location. If one is concerned with turbu- 
lence models to be used for simulator studies of STOL airplanes 
in approaches to airport, during which altitude varies, the 
lack of homogeneity becomes important since the statistical 
properties vary with altitude. 

Ergodicity . - Since, in turbulence measurements, it is 
virtually impossible to obtain ensembles of turbulence data 
from atmospheric measurements, it is necessary to use time 
averages to get statistical information. If such a time 
average gives the same statistical properties as the ensemble 
average the process is called ergodic. Note that ergodicity 
implies stationarity. 

In the various methods of turbulence simulation it is not 
always necessary to assume either stationarity or ergodicity 
since it may be possible, if desired, to vary the statistical 
properties with time. However, in the reduction of experi- 
mental data, whether turbulence, flight, or flight simulation 
data, only a few relatively short time histories are generally 
available, and as a result one usually assumes ergodicity. 

Note that all the quantities defined below employ time 
averages and therefore assume ergodicity. 

Average or Mean Value . - The mean value of an ergodic 
random process is 


Mean of u 


u 


= ilim 

T^OO 


i /■ 


2T 


u(t) dt 


-T 


( 1 ) 


For practical evaluation of the mean value the limit process 
is not required and u can be approximated by the equation 
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(2) 


- 1 

u ~ ^ /q u(t) dt 


where T is large. 

This approximate representation is especially useful for 
processes such as turbulence which are probably not stationary 
over very large time periods. However, the time interval T 
must be large enough so the averages approach the asymptotic 
value one would obtain for a stationary process, but small 
enough to be assured that the process has approximately con- 
stant statistical properties during this time interval. 

Variance.- The variance of u is defined as 


Variance of u = a ^ = iim ^ / [u(t)-u]^dt (3) 

ip^OO ^rji 


As before for practical purposes the variance can be approxi- 
mated by 


T 

""u^ ~ k ^0 dt (4) 


for sufficiently large T . 

Standard Deviation, Root Mean Square .- The standard 
deviation is defined as the square root of the variance, or 


Standard deviation = a = /a , (5) 

u u ' ' 


which is the root mean square of the departure from the mean 
value. In this report it is assximed that u is equal to 
zero, or that signals are measured from the mean value. Thus 
the standard deviation becomes the root mean square of the 
signal itself. 
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4“ 

Normalized Central Moment.- The N normalized central 

moment/ or more simply the N moment/ of a random process 
u(t) is 


M 


n 


ilim 

•p->-oo 


2T 


T 

f 

-T 


... “1 

!iM^ dt 

a 

U u J 


n 


1 , 2 


3 . . . 


(6) 


or approximately 


M ~ 
n 


1 

- f 
T 0 


u (t) -u 
a 

u 


n 


dt / n = 1 / 2 / 3 


(7) 


The moment is called central because the mean value u is 
substracted from u(t) / and is nondimensionalized with res- 
pect to u by dividing by the standard deviation . 

Inspection shows that M.. is thus zero, and becomes 

a 2 

< 5 ^> = 1 • 

u 

Cumulative Probability Distribution .- The cumulative pro- 
bability distribution is defined as 


3P^(x) = Probability that u < x 


( 8 ) 


Two important properties of the cumulative probability are 
that IP (-“) = 0 / which implies that no value of u will 

be less than -" , and (+“) = 1 / or that all values of 

u will be i +“ . If a random function is distributed 
symmetrically about u = 0 / that iS/ there is just as much 
probability of having a negative signal of a certain size as 
a positive signal/ (0) will be 0.5 , and the cximulative 

probability distribution will look about as shown below. 
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1 . 0 1 


Cumulative 

Probability 



Probability Density Distribution .- Since for random pro- 
cesses it is difficult to conceive of the probability of an 
exact event occurring i such as having a g'ust of exactly 
10 m/sec (such a probability would be almost zero ) , the pro- 
bability density distribution is defined as. 


Probability density = P^(x)dx = Probability that x < u < x + dx 

(9) 

The probability density is related to the cumulative probability 
by differentiation. 



will be greater than <*> are zero, and 
00 

/ P (x)dx = 1 , or the probability that u 

—00 

will lie between ±«> is 1 . 
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Gaussian Probability Density Distribution .- The most 
commonly known probability density distribution is the Gaussian 
or normal distribution 


P„ (X) = — ^ exp [-1/2 (^) ] 


( 11 ) 


u 


a /Tn 
u 


u 


When u is assumed to be zero, as it is in this report ^ 


P, , (X) = — ^ exp [-1/2 (^) ] 

^ a /2 tt °u 

u 


(12) 


Examples of the Gaussian distribution for two values of 

a are shown below, where u = 0 . 
u 


0.1 

P(x) 





Modified Bessel Probability Density Distribution ♦- A 
probability density distribution which will be referred to 
frequently in the following report is characterized by the 
modified Bessel function of the second kind and order zero, 


K. 




u 




( 13 ) 
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A comparison of the modified Bessel function probability 
density and the Gaussian probability density for a = l m /sec 
is shown below ^ 



X 


An examination of these probability density distributions 
shows that; 

1. ) The Bessel function shows greater probability for 
small gusts (0.5n/sec or less) than does the Gaussian. 

2. ) Between gusts of about 0.5m/sec and 2 m/sec the 
Bessel function indicates less probability of occurrence 
than does the Gaussian. 

3. ) The Gaussian shows almost no probability of reaching 
gust velocities of over 3 m/sec, while the Bessel func- 
tion shows much greater probability of reaching such 
higher gust velocities. 

Higher Moments .- As mentioned on page 15, the moment 
of a process in whxch u = 0 is 


M ~ ^ / 
n To 


(Ultl, 

u 


n 


dt 


n = 1 


(14) 


As indicated the 1st moment is zero, and the second moment 
becomes 1 . The higher moments are important in the des- 
cription of the probability density distribution. The third 
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moment is called the "skewness" moment, and if a process is not 
symmetrical about the mean value the third moment will be non- 
zero as shown below. 

P(x) I 



Since in this report all processes are assumed to be symmetric 
about, the average, the third and higher odd moments will be zero 

The fourth moment. 



/ U(t) . 

' a ’ 
u 


( 15 ) 


is called the kurtosis or flatness factor, and probability 
densities with small and large kurtosis could look as below, 
for the Scime value of o . 




Comparison of these curves shows that the higher (even) moments 
are very important in defining the shape of the probability 
density, and that probability densities with larger kurtosis 
have higher probability of obtaining large values of the gust 
velocity. It will be noted that this same characteristic has 
been pointed out on page 18 in regard to the Bessel function, 
which has higher probability of having large gusts than does 
the Gaussian probability density. 

Cross Correlation Function .- The cross correlation function 
of two random processes u (t) ^and w(t) is 


1 

C (t) = lim f u(t)w(t+x) dt 

-T 


(16) 


Correlations are the measure of the predictability of a signal 
at some future time (t + x) based on the knowledge of a signal 
at the time t . Since the modeling of a random process such 
as turbulence consists of developing techniques for predicting 
the time-wise behavior of that process it can be seen that 
correct duplication of the correlation functions is very impor- 
tant, since these are measures of predictability. 

Autocorrelation Function .- The autocorrelation function of 
a random process is a special case of the cross correlation 
function defined above in which w(t) = u(t) , so 


1 

C„„(x) = lim ^ / u(t)u(t + x) dt (17) 

-T 


Integral Scale Length .- A statistical parameter of special 
importance in atmospheric turbulence studies is the integral 
scale length, also called simply the scale length. 


U 00 

L = %r f C (x) dx 

u „ 2 ^ uu' ' 

O “"00 
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(18) 


where u^ is the reference steady state flight speed of the 
aircraft^f lying through the turbulence. 
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In random process theory an integral scale time is 
defined as 


J 


00 


/ 


C (t) 
uu 



dx 


(19) 


As shown below, this value of J is a measure of the strength 
of the correlation, and gives a rough idea of the time interval 
during which the process is correlated with itself. 



In work concerning the flight of an aircraft through turbulence 
it is more convenient to turn this time into a distance by 
multiplying by the steady flight speed, or 

L = u j" , (20) 

u 0 

so the scale length becomes an approximate measure of the dis- 
tance an aircraft flies through turbulence which is correlated 
with itself. 

Cross Spectral Density .- The cross spectrum of two random 
processes u(t) and w(t) is defined as the Fourier transform 
of their cross correlation. 


00 

(f) = / C,„(x) exp{- i2Trfx) dx 

uw uw 


( 21 ) 


where f = frequency in hertz . 


21 



The Fourier Transform of an aperiodic function such as 
can be interpreted as a continuous spectrum of Fourier 

series coefficients for that function. Thus <l> is a 

harmonic representation of C . 

uw 

The cross spectral density of two random processes is 
generally a complex function of frequency. The phase angle at 
any frequency represents the average phase difference of the 
Fourier components of the two processes at that frequency. 

It is also noted that in this report, for mathematical 
convenience. is a function of both positive and nega- 

tive f , although negative f has no physical significance. 

Power Spectral Density .- The power spectral density of 
a random process is the Fourier transform of its autocorrela- 
tion function, or 


<j> (f) = f C (T) exp(- i2TTfT) dT (22) 

U UU 


Unlike is always real since there is no phase 

information contained in an autocorrelation function. 

<f> can be interpreted physically as the average contri- 
“ 2 

bution to the variance from the frequency component f . 

The variance of u(t) can be obtained from by integration. 
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White Noise .- White noise is a random process for which 
the power spectral density is a constant, independent of 
frequency. That is 


= K (24) 


where K is a constant for all f . True white noise cannot 
be realized physically since its variance 
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= / (p (f) df 
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would be infinite. However, random processes for which the 
power spectrum is constant over a broad frequency range can 
be generated. This is called band limited white noise. When 
band limited white noise is used as the input of a system 
which responds only to frequencies in the constant power band 
of the noise, then all of the theoretical results obtained by 
assuming true white noise are valid. 


Description of Atmospheric Turbulence 

Before specifying the statistical properties of a realistic 
turbulence simulation it is necessary to review the characteris- 
tics of actual turbulence. In this report the statistical pro- 
perties which are judged to be of most importance are the pro- 
bability distribution, power spectral density, cross spectral 
density, and patchy characteristic of the gust field. At this 
time only the gust components occurring at one point in space 
will be considered. Later in this section, when discussing the 
method of turbulence model input to the simulator, the proper- 
ties of real turbulence described here will be expanded to 
include the spatial characteristics of the gust field. 

Probability Distribution .- The probability distribution of 
a random function provides information concerning the range of 
values assumed by that function and the frequency with which 
they occur. There are several ways to present such information. 
The methods described here are the cumulative probability dis- 
tribution (eq. (8) ), the probability density distribution 
(eq. (9) ), and the normalized central moments (eq. (14) ). 

The cumulative probability distribution and the probability 
density of a random process are related by differentiation as 
shown in equation (10) , and therefore contain identical informa- 
tion. However, this information is presented by the two distri- 
butions so as to make apparent different aspects of the random 
process. The probability density function, usually plotted on 
linear scales, is useful for comparing differences in the pro- 
bability of small values of a random process. The cumulative 
distribution on the other hand, especially when plotted in pro- 
bability coordinates (see fig. 1) shows the behavior of extreme 
values of the process. 

A third way to present the same sort of information is 
through the use of normalized central moments. Though these 
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moinents do not contain the amount of information present in the 
cumulative probability or probability density functions, they 
are much easier to estimate and do give some indication of the 
occurrence of large values of the random process. 

As there is little experimental data available which distin~ 
guishes between probability distributions of the different gust 
components, no distinction will be made in this report. That is, 
all components of the turbulence model will be assumed to have 
the same probability distribution. 

Cumulative Probability Distribution .- Figure 1 presents 
data from reference 5 showing typical cumulative probability^ 
distributions of atmospheric turbulence. The straight line in 
that figure represents the result expected from a Gaussian 
process. The departure of the turbulence data from the Gaussian 
curve clearly indicates the increased probability of large gusts. 
Data are shown for two different altitudes, and similar tenden- 
cies are seen in both cases. An analysis of unpublished data 
collected by the University of Washington Department of Atmos- 
pheric Sciences from a four meter tower during times of strong 
winds, and therefore nearly neutral stability conditions, indi- 
cates that data of the form shown in figure 1 is indeed repre- 
sentative of low altitude turbulence. 

Probability Density Distribution.- A typical probability 
density (eq. (9) j obtained from low altitude turbulence is 
compared with the Gaussian distribution in figure 2. The 
increased probability of large and small gusts in real turbu- 
lence is apparent. 

Normalized Central Moments . — The table below compares the 
fourth and sixth normalized central moments (eq. (14) ) measured 
in low altitude turbulence with those expected from a Gaussian 
process. These data were obtained from LO-LOCAT Tests, Phase 3, 
Test 1141. The larger values measured in real turbulence indi- 
cate the presence of larger gust velocities than would be 
expected from a Gaussian process. 

TABLE 1.- COMPARISON OF FOURTH AND SIXTH NORMALIZED CENTRAL 

MOMENTS OF TURBULENCE WITH THOSE OF A GAUSSIAN PROCESS 




«6 

Atmospheric Turbulence* 

3.5 

21.7 

Gaussian Process 

3.0 

15.0 
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*Data from 76.2 m altitude, neutral stability, 
strong winds 




It should be noted here that no attention has been given to the 
odd moments of turbulence such as M, or M . Experimental 
measurements of laboratory generated turbulence (ref. 6) indi- 
cate that these odd moments are not zero. However, there does 
jjot appear to be any detailed data on these moments measured in 
atmospheric turbulence. In any case there does not appear to 
be any convenient way to introduce these odd moments into a 
turbulence simulation and they have been neglected in this report. 

Power Spectral Density .- The power spectral density of a 
random process (eg. (22) ) provides information on the average 
contribution to the process from the frequency components which 
make it up. It is therefore a most important statistical des- 
cription. Several possible algebraic forms have been suggested 
for the power spectral density of the three turbulence gust 
components. A comparison of the most popular forms presented 
in reference 7 indicates that the most accurate form is that 
suggested by Theodore von Karman (ref. 8) to represent isotropic 
turbulence 
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These shapes are not convenient for turbulence model work because 
they cannot be exactly matched using linear filters. (See dis- 
cussion of linearly filtered Gaussian white noise turbulence 
simulation which follows in this section.) This is because of 
the non-integer powers appearing in the denominators. The error 
resulting from the use of linear filters to produce these spec- 
tra can be reduced as far as desired by adding poles and zeros 
to the filter transfer function, but this adds to the complexity 
of the system and is generally felt to be not worth the trouble. 
This problem is usually overcome by assuming the spectra sug- 
gested by H. L. Dryden (ref. 9). 
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These spectral shapes can be exactly matched using linear filters. 
Since the difference between the Dryden shapes and those proposed 
by von Karman are quite small in the frequency range of interest 
in most flight simulation studies, there is no reason to believe 
that these spectra are not sufficiently realistic. 


In this report a slightly modified form of the Dryden spec- 
tra will be used. The changes introduced are: (1) replacement 

of the single scale length L with independent values for each 
gust component. This change is consistent with much experimental 
data (ref. 10 for example) which shows that the scale length 
does vary among the components of atmospheric turbulence; (2) 
replacement of the lateral power spectrum 4> with a form 

similar to that used for the longitudinal spectrum 4>^ • This 

change is suggested by low altitude data (ref. 10) showing the 

lateral spectrum to be more like 4> than like ^ . It is not 

u w 

a strict requirement of the turbulence simulation schemes to be 
described presently, but does simplify them slightly. 


The gust component power spectral densities to be used in 
this report are: 
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These spectra are presented non-dimensionally in figure 3. 

It should be noted that these spectra 
for anv heading of the simulated vehicle relative to the mean 
wind! in SiSLl the scale length and standard deviation of 
each gust component are expected to be functions of the 
angle! altituL, atmospheric stability 

tive to the surface, though explicit relationships are not 
needed for the work reported here. Many reports describing 
the dependence of L and a upon the above 

tions with the exception of heading angle, can be found in the 
ci?!eAt Utera!ure if explicit functional relationships are 
?equi!ed. (See ref. 11 for example.) The dependence upon 
heading angle is generally unknown. 

cross Spectral Densities .- The cross spectral density 
(eq (21) ) of two random processes provides information on the 
pha^e relationships of their frequency conponents At the pre- 
Lnt time little is known about the cross 

of the u , V , and w gust components, particularly the varia 
tion of these cross spectra with altitude, r „®erest 

tive to the mean wind direction, and other fSrms 

in fliqht simulator work. Because convenient algebraic forms 
for these spectra are not currently available, and because 
their inclusion in a turbulence simulation , , 

plicates the simulation technique, cross spectra are ^ ^ 

nealected in flight simulations. However, one would expect 
“’llelrdSss spectra in any turbulent boundary layer |l°w- 
and therefore in most low altitude turbulence. Ind«d Elderkin 
fref 10) finds a consistent low frequency cospectrum |reai 
pllt orthe Somplex cross spectrum) relating the longitudinal 
and vertical gust components measured from a tower. 

In aeneral it would seem prudent to make some provision 
for loS f?equencrcross spectra in a turbulence simulation even 
no specif iS form is known. In this report a low frequency 
cross spectrum relating the longitudinal and lateral gust com 
pSnSntJ'^wIn“e assume! The form of this spectrum is determined 
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primarily by the technique in which the artificial turbulence 
time histories are to be generated, and therefore involves 
parameters which are dependent upon the simulation method itself. 
Details of the cross correlation method used in this report and 
the physical significance of the parameters will be found in 
Appendix 1. The cross spectral density assumed in this report is 
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where A , B , and R are arbitrary parameters satisfying the 
conditions 
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Patchy Nature .- It is known that turbulence has a patchy 
structure. That is, atmospheric turbulence seems to occur in 
patches of relatively intense motion separated by areas of rela- 
tive calm. It is not easy to formulate a model of this phenom- 
enon in terms applicable to flight simulator work. Nor is it 
easy to see a patchy structure in turbulence velocity time his- 
tories, although the derivative of a gust velocity time history 
as shown in figure 4 does reveal distinct bursts or patches of 
activity. (Reference LO-LOCAT Phase 3, Test 1145.) Since there 
are no available models of this patchy characteristic, figure 
4 will be used as a standard of comparison in this report. 
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Presently Used Simulation Techniques 

The preceding pages of this section have discussed the 
statistical properties of atmospheric turbulence which are to 
be modeled by a realistic turbulence simulation. In this sec- 
tion several presently used simulation techniques will be des- 
cribed from the standpoint of their statistical realism and 
suitability for use in flight simulators. 

Flight Recordings of Atmospheric Turbulence .- Perhaps the 
most obvious method of producing a realistic simulation of 
atmospheric turbulence is to use recorded time histories obtained 
under the same flight conditions as those being simulated. There 
can be no question as to whether or not these time histories are 
an accurate representation of atmospheric turbulence. However, 
this method has several serious drawbacks. First, the time 
history must be obtained under the scime conditions as are to be 
simulated. This means that variables such as altitude and true 
airspeed must be matched, as well as meterological variables. 

This implies that a flight simulator facility which is to 
vide realistic simulations of a wide variety of flight conditions 
must have a complete library of time histories available. A 
second problem associated with the use of recorded time histories 
is that the length of each time history is fixed. Longer run 
times cannot be accommodated without an interruption of the 
simulated turbulence. One further problem with recorded time 
histories is that repeated use of the same time histories may 
permit the pilot to become familiar with some of their charac- 
teristics, such as the occurrence of large gusts. This would, 
of course, negate the principle purpose of the artificial turbu- 
lence, which is to provide an unpredictable external disturbance 
to the flight simulation. 

The recorded time history turbulence model does have advan- 
tages in certain circiimstances . For example, it can be used to 
simulate statistically rare occurrences such as encounters with 
isolated large gusts. Since this event is very infrequent, a 
turbulence simulation employing the statistics of typical turbu- 
lence, described earlier in this section, would produce such an 
encounter only very infrequently. This is, of course, an 
unacceptable situation for a practical flight simulator program. 

In summary, the recorded time history turbulence simulation 
does have application to the simulation of situations which are 
not adequately described by the statistics of typical turbulence 
discussed earlier in this section. However, it lacks the versa- 
tility to simulate a wide variety of flight conditions, is 
limited in the length of its time histories, and may permit the 
pilot to become familiar with time histories which are used 
repeatedly. It is, therefore, not suitable for the simulation 
of typical turbulence. 
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Sum of Sine Waves.- Another technique which has been used 
to simulate atmospheric turbulence is the sum of sine waves 
model. This technique involves summing ten or twenty sinusoidal 
waves of different frequencies and magnitudes. The resultant 
time history is then used to represent a time history of turbu- 
lence. This method does provide a random appearing, non-repeat- 
history, but cannot be said to represent atmospheric 
turbulence, which contains an infinite number of frequency com- 
ponents. If the number of sinusoids is increased in an attempt 
^ smoother power spectrum, the resultant time history 
Will become Gaussian by virtue of the central limit theorem. 

As discussed previously in this section, atmospheric turbulence 
IS not Gaussian. Hence the sum of sine waves turbulence simula- 
tion IS not statistically representative of atmospheric turbulence. 

^ Lij^early Filtered Gaussian White Noise .- The most commonly 
used turbulence simulation is the linearly filtered Gaussian 
white noise technique. Each component of this simulation is pro- 
duced by passing Gaussian white noise through a linear filter as 
shown in figure 5. As noted in that figure, the power spectral 
density of the simulated time history is determined by the trans- 
fer function of the filter. Any power spectral shape can be 
approximated as closely as desired by using a sufficient number 
of poles and zeros in the transfer function. 


This turbulence simulation is remarkably easy to implement 
and can be generalized (see for example reference 12) to produce 
several simultaneous random time histories with specified power 
spectral densities and cross spectral densities. Furthermore 
these spectra and cross spectra are continuously variable by 
changing the parameters of the filter transfer functions. 

Despite these desirable characteristics, this type of turbu- 
lence simulation also has its deficiencies. Since a Gaussian 
process r^ains Gaussian when passed through a linear filter, 
the artificial turbulence time histories produced by this method 
will always be Gaussian. As already noted, atmospheric turbu- 
lence IS not Gaussian. 

It is also easily shown that the Gaussian white noise turbu- 
lence model does not reproduce the desired patchy nature of 
turbulence. Figure 6 compares the derivative of an artificial- 
turbulence time history produced by the Gaussian model with the 
typical example of atmospheric turbulence from figure 4. Care 
was taken to match the power spectral densities of the two 
curves, so both time histories have the same frequency content. 
Note that the Gaussian time history completely lacks the patchy 
nature of real turbulence. 
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In addition to the criticisms based on the poor statistical 
characteristics of the Gaussian turbulence simulation, there 
are also indications that this method produces time histories 
which do not "feel" like turbulence when used in a flight simula- 
tor. Pilots complain that compared to real turbulence, the 
artificial turbulence is too regular and does not require correc- 
tive control motions, because a disturbance in one direction is 
soon followed by a compensating disturbance in the opposite 
direction (ref. 13) . It must be concluded therefore that the 
artificial turbulence generated by the linearly filtered Gaussian 
white noise turbulence simulation is not sufficiently realistic 
for flight simulator applications. 

In summary, the Gaussian turbulence simulation has many 
convenient features. It is easily generated and can be gener- 
alized to produce a number of simultaneous time histories with 
specified power spectra and cross spectra. However, the method 
does not model the non-Gaussian nature of atmospheric turbulence, 
nor does it model the patchy character of real turbulence. It 
is therefore unsuitable for realistic flight simulator applications. 

Method of Orthogonal Functions.- A recent report by Dutton 
and Deaven (ref. 14) describes a method by which turbulence time 
histories are decomposed into eigenfunctions of the covariance 
matrix. Simulated turbulence time histories are then generated 
by recombining these eigenfunctions. Only preliminary results 
are presented, however it appears that even though the - 5/3 
slope of the power spectral density is correctly modeled, the 
procedure as now proposed does not correctly model the non- 
Gaussian nature of turbulence. The method also shares some pro- 
blems of the recorded time history simulation in that the time 
history length is fixed, and maneuvers of the simulated vehicle 
cannot be correctly simulated unless they are known in detail 
before the simulator experiment, and the appropriate data have 
been collected. The primary advantage of this method is that 
the probability of occurrence of any particular time history 
can be computed a priori and therefore the degree to which the 
results are representative is known. In summary, the method of 
orthogonal functions almost certainly has flight simulator 
applications in the area of rare event simulation, such as the 
isolated large gust encounter described previously in this sec- 
tion, but does not provide the versatility being sought in 
this report. 


A Proposed Non-Gaussian Turbulence Simulation 

The preceding parts of this section have discussed the pro- 
perties of atmospheric turbulence and several simulation techni- 
ques. Of the simulations described, the linearly filtered Gaussian 
vhite noise method is the most easily produced and is the method 
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most readily adaptible to simulating complex situations involving 
time-varying power spectra and cross spectra. Unfortunately the 
artificial turbulence which it produces is always Gaussian and 
is therefore unrealistic. It seems reasonable to suppose that a 
turbulence simulation similar to the Gaussian method, but with 
a non-Gaussian probability distribution, will be more realistic. 

A simulation of this type will now be proposed. 

Consider the patchy nature of turbulence. The existence of 
patches suggests that atmospheric turbulence might be simulated 
by a product of Gaussian processes, the first representing a 
constant-intensity turbulence time history and the other repre- 
senting the time varying standard deviation. This idea has 
already been pursued in reference 4, and a non-Gaussian turbulence 
simulation has been proposed. Furthermore, reference 4 shows 
that this new simulation has the same versatility as the Gaussian 
model. That is, it is possible to generate the reguired power 
spectra and cross spectra, with time-varying characteristics to 
represent vehicle maneuvers just as in the case of the Gaussian 
turbulence simulation. However, the Gaussian product simulation 
is not as realistic as might be hoped, as will be discussed later. 

The probability density resulting from a product of Gaussian 
processes is the modified Bessel function distribution defined 
previously in this section (eq. (13) ) . This distribution is 
compared with the normal distribution and a typical result for 
atmospheric turbulence in figure 7. Note that the probability 
density of real turbulence lies between the Gaussian and Modified 
Bessel distributions. This suggests that a more realistic pro- 
bability density would be obtained by adding the two. 

Typical patchy characteristics of the Modified Bessel pro- 
cess are compared to those of the Gaussian model and an actual 
turbulence time history in figure 8. All three functions shown 
in figure 8 have the same power spectral density and therefore 
the s^e frequency content. However, the Modified Bessel pro- 
cess is clearly far more patchy than atmospheric turbulence. 

Table 2 below presents the fourth and sixth normalized 
central moments of Gaussian, real, and Modified Bessel processes, 
and shows that the Modified Bessel greatly overestimates the 
higher moments, while the Gaussian underestimates them. This is 
further evidence that a summation of Gaussian and Modified Bessel 
processes might produce a more realistic result than either of 
these above. 
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TABLE 2.- COMPARISON OF FOURTH AND SIXTH NORMALIZED CENTRAL 
MOMENTS OF TURBULENCE WITH THOSE OF GAUSSIAN AND MODIFIED 
BESSEL PROCESSES 



^4 

^6 

Gaussian 

3.0 

15.0 

Real Turbulence 

3.5 

21.7 

Modified Bessel 

9.0 

225.0 


*Data from 76.2 m , neutral stability, 
strong winds. Ref. LOLOCAT, Phase 3, 
Test 1141 


The large moments of the Modified Bessel process indicate 
that it has much larger peak values than are found in typically 
measured low altitude turbulence. The Modified Bessel model 
may therefore be suitable for simulating some extremely violent 
flight environment, but it not very suitable for representing 
typical turbulence. 

The non-Gaussian turbulence simulation proposed ® . 

report, which consists of the addition of Gaussian and Modified 
Bessel processes, is shown schematically in figure 9. The 
detailed analysis of this system and its application to a spe- 
cific turbulence simulation is both complicated and lengthy. It 
has therefore been placed in Appendix A. The following is a 
summary of some of the proposed simulations important features. 

In figure 9 a(t) , b(t) and d(t) are Gaussian processes Pto- 
dSced just as in the case of the linearly filtered Gaussian white 
noise turbulence simulation. Therefore a(t) , 

are independent Gaussian processes. The process c(t) is je^er 
ated by the multiplication of the Gaussian processes a(t) an 
b(t) . it is therefore a Modified Bessel process. The system 
output, u(t) , is the sum of the Modified Bessel process c (t) 
and the Gaussian process d(t) . 

The probability distribution and patchy character of u(t) 
are determined by the standard deviation ratio 


a 


R = 


(39) 


If R is very small then a is much smaller than , and the 
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process c(t) can be neglected in comparison to d(t) . The 
system of figure 9 is then equivalent to the Gaussian simulation 
shown in figure 5, and the system output is therefore Gaussian. 

The characteristics of the output are then the same as those of 
the Gaussian process in figure 8 and Table 2. On the other hand, 
if R IS very large then is much larger than cr^ , The 

Gaussian process d(t) can be neglected and the system output 
IS the Modified Bessel process. The characteristics of the 
output are then the same as those of the Modified Bessel process 
in figure 8 and Table 2. 

It can be seen that the standard deviation ratio R is a 
powerful control over the probability distribution, and that a 
close approximation of almost any reasonable experimentally deter- 
mined distribution is available. It has been found that, in 
general, a good match with measured probability data is obtained 
when R is of the order of unity, which results in a random pro- 
cess with properties lying between those of the Gaussian and 
Modified Bessel processes. As mentioned above, the theoretical 
analysis of the proposed non-Gaussian turbulence simulation and 
Its application to the generation of specific turbulence models 
IS quite lengthy and has been placed in Appendix A . 

Some theoretical results showing the dependence of the system 
oi^tput characteristics upon the ratio R are presented in the 
following mentioned figures in order to show the range of varia- 
tion possible with the proposed model. Figure 10 presents the 
range of cumulative probability distribution which can be obtained. 
Note upon comparison with figure 1 that the probability of large 
gusts can be varied over a range which includes the distribution 
of real turbulence. Figure 11 presents the probability density 
distributions corresponding to the cumulative distributions of 
figure 10. The experimental data from figure 2 have been plotted 
in figure 11, and it is evident that the proposed model with R 
approximately unity compares very favorably with these data. 

Figure 12 shows the dependence of M^ and Mg upon R . Figure 

13 presents the time derivative of the system output for several 
values of R . All four time histories have exactly the same 
frequency content. Note that the patchy characteristics of tur- 
bulence as typified by figure 4 can be matched quite closely. 


10 through 13 show that the various properties of 
can be matched quite well if R is chosen properly 
in each case. Of course these properties are not independent of 
each other. They all depend upon R and only one value of R 
can be chosen for any given turbulence simulation. It is there- 
fore necessary to choose a value which is in some sense best. 

A comparison of the above mentioned figures with the data pre- 
sented for real turbulence in figures 1, 2, and 4 and Table 1 
indicates that an R value of slightly less than unity appears 
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to give very good results. Of course the final test of realism 
cL only coL from a practical flight simulator exper iment . As 
yet no Attention has been given to the explicit 
the artificial turbulence time histones should be introduce 
into the simulator so as to provide a realistic representation 
of flight through turbulent air. This is the next sub] 
this report. 


Simulation of Spatial Distribution Effects 

The simplest way to handle gust inputs is to assume the 
gusts are uniform everywhere over the aircraft and are equal to 
the gusts measured at the aircraft center of gravity • 

These gusts are then fed into the airplane equations of motion 
through stability derivatives such as ^ ' • • f 


^3 


m 


a 


m 


= C a 

III 

a 


+ c + 

m u^ 
a 0 


This method is called the one-point model, and is in fact quite 
accurate for moderate and high altitudes, and away 
shear surfaces, where the ratio L/b of scale-length-to-span 
is vLriLge (say > 20) . (L can be loosely thought of as 
being approximately the size of the average turbulent eddy. For 
wind^driven turbulLce in the lower part of the atmosphere, 

L?h where h is the altitude (ref. 15). For higher alti- 
tudes, 'above 533 m (1750 ft.), L :: 533 m) . For large values of 
L% the eddies are very large compared to the wing span, and 
the gusts are very nearly uniform over the whole airplane. As 
Vb drops below 20, gradients in the gusts become 
across the wing span and from the c.g. to the tail, and these 

gradients cause forces and moments ITlne 

dered bv the one-point model. For example, an up-gust at one 
w!ng-tip ani a down-gust at the other causes a 
on ?he Lrcraft. Also, an up-gust at the ^ 
be different from the up-gust at the horizontal tail. Th 
would result in a nose-up or nose-down pitching moment not 
considered in the one-point model. The correct simulation of 
these effects is very important for a STOL airplane at lo 
S?i?ude In this report, the De Havilland Twin Otter (with 
a wing span of 19.8 m (65 ft.) ) was sunulated at an altitude 
of 76.3 m (250 ft.). For this case L/b ~ 4 , 

Several theories have been proposed to account for the 
effects of the spatial distribution of gusts over an aircra . 
Itkin treL 16 Tn) expands the gusts in a Taylor series about 
the S.g! using a second-order series, as high as one can go 
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before the method becomes unwieldy, results in the gust model 
being accurate for L/b down to 6 (too high for the case 
considered in this report). Skelton (ref. 12) proposes a 
3-point gust model, with the gusts considered as acting at the 
tail and mid-way along each wing. His method requires eight 
separate gust inputs, all properly cross-correlated, and still 
only depicts the gusts as being linearly distributed across 
the span. This method would be accurate for L/b down to 
only about 12. Eggleston and Diederich (ref. 18) use a modified 
strip theory to calculate the power spectrum of the rolling and 
yawing moments on an aircraft wing in turbulence. Their method 
is good for any ratio of L/b , although they account only 
approximately for the unsteady aerodynaunic effects. Franklin 
(ref. 19) makes use of their results in his flight-simulator 
investigation into lateral-directional flying qualities. 


The method of handling the gusts in the present simulation 
generally follows Franklin's ideas, with the exception of three 
differences : 


1. ) The method of handling tail lags in that Franklin uses 
a linear-filter approximation to the Laplace transform of 

a time delay, e , whereas the present simulation uti- 
lizes a separate gust input at the tail; 

2. ) This report includes the effects of distributed u 

g 

gusts across the wing span (a necessary inclusion since a 
slow-flying STOL aircraft is flying at a large angle of 
attack) ; and 


3.) The inclusion by Franklin of more unsteady-aerodynamic 
terms, by multiplying the , and w^ input power 

spectra by a Sears function (ref. 20) . 

^ In the present simulation, the gust components were handled 
in the following way; 


a.) 


Side gust component, 


V 

g 


(1) Gust at airplane c.g.; v (t) 

^c 

This appears in the sideslip angle 
by the stability derivative of the 
in the equation for yawing moment. 


3 and is multiplied 
tailless airplane, e.g 


/ 


V 

c = (c - c ) -3£ + 

^ ^3 ^t3 '"rw 


where 


36 



C = the yawing-moment-due-to- sides lip stability 

^3 derivative for the complete airplane, 

C = the yawing-moment-due-to- sideslip stability 

*^t3 derivative due to the tail only, 

V /V = the effective sideslip angle measured at the c.g. 
gc RW 

(2) Gust at tail (gust at c.g. delayed): 


' ''g ■ ir > • 

y t o 


This accounts for the finite time it takes for a gust to 
travel the distance between the wing and the tail. 

(The frozen— field hypothesis is assumed for the gust field, 
i.e., the gust field is assumed fixed, unchanging in space, 
while the airplane flies through.) The time-lag ~ 

.212 seconds was used here, v (t) appears as a sideslip 

^t 

angle at the tail, and influences the equations through 
the tail contribution to the sideslip stability derivatives 
such as . The complete yawing moment coefficient 

due to gust^^then becomes 


n 


= (C 


n, 


- c 


n 


t6 


V. 


_ 2 £ 


+ 


RW 


n 


tB 


'^RW 


The tail-only stability derivatives were estimated using 
methods in references 21 and 22. The complete equation 
for the yawing moment is shown in Appendix B. 

b.) Horizontal and vertical gust components, Ug and w^ 

(1) Gusts at airplane c.g.: u (t) and w (t) . 

^c ^c 

These influence the equations of motion through the tail- 
off stability derivatives with the w gust appearing as 

^c 

an effective a (and corresponding a term) . 
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(2) Gusts at tail (gusts at c.g. delayed) : 




(41) 




and 


(42) 


These appear through the tail contribution to the stability 
derivatives, i.e., C and c 

"*ta "‘tu 

(3) Rolling moments due to spatial effects on wing- 
(t) and C. (t) 

ug ^wg 

These appear as additive terms in the rolling moment equa- 
txon, and are generated by taking two independent random- 
noise signals and passing them through appropriate filters 
to shape the power spectra to match those of Eggleston and 
Diederich (ref. 18), including the unsteady-aerodynamic 
effects. (See page 39.) 

(^) Yawing moments due to spatial effects on wings 
(t) and C (t) . 

ug "wg 

These are again additive terms in the yawing-moment equa- 
tion. Reference 18 shows that these terms are simply a 
constant multiplied by the rolling -moment spatial effects; 


ug ug 

(t) = ^ (t) (44) 

wg Jlp wg 

In summary, the turbulence model consisted of the following; 
gust inputs at the c.g. and acting on the tailless air- 
plane; the same gusts delayed by ^^/u^ and acting on the tail; 

and random rolling and yawing moments to account for the distri- 
bution of gusts over the wing span. Required for the simulation 
are three time delays (done digitally here) plus time histories 
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of the five inputs u (t) , v (t) # (t) / C« (t) , and 

Q q q JC _ 

ug 

(t) . In this simulation, these inputs were computed before- 
wg 

hand, recorded on FM tape, and input to the simulator as analog 
signals . 


Power Spectra of and 

ug wg 

The power spectra of the spatial effects (t) and 

ug , 

(t) use the results of reference 18. In that report, 
wg 

Eggleston and Diederich considered a number of different wings 
in turbulent flow. These wings were such that they had differ- 
ently shaped steady-state lift distributions for the steady- 
flow case. For the turbulent case, the authors used quasi- 
steady aerodynamics (in which the lift distribution is assumed 
to adjust instantaneously to changes in the local angle of 
attack) to derive power spectra for the rolling and yawing 
moments on the wing. Surprisingly, the resulting non-dimensional 
spectra turned out to be practically independent of the 
lift distribution assumed. Eggleston and Diederich suggested 
taking into account the unsteady-lift effects (i.e., the fact 
that the lift does not respond instantly to a change in local 
a ) by simply multiplying the power spectral density of the 

2 

i,\n 

rolling moment due to each gust component by I S, 


(— ) 
'k '2uJ 


where 


is the Sears function (ref. 20) . That was what was 


done here. The present simulation takes the power spectra from 
figures 7b and 9b of reference 18 for the elliptic span-loading 
case and b/L = 0.25 . Since the span b of the Twin Otter 
is 19.8 m (65 ft.), this corresponds to a turbulence scale- 
length of 79.4 m (260 ft.), which corresponds roughly to the 
altitude of 76.3 m (250 ft.) used in the simulation. The result- 
ing power spectra of rolling moments, and 4>p (in 

ug wg 

the present notation) are shown in figures 14 and 15 respec- 
tively. Unsteady-lift effects are included by using the Sears 
function as outlined above. Note that in converting from 
Eggleston and Diederich' s curves to the present simulation, 
the value = - .337 was used for the zero-lift angle of 

attack of the wing. Note also that in reference 18 the power 
spectrum is defined as 1 /tt times the way it is defined in 
this report. 
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For convenience in generating the time histories of 

(t) and Cj^ (t) , Eggleston and Diederich's power spectra 
ug wg c- 

(with unsteady-lift effects included) were approximated by a 
Ixnear-f liter representation, i.e. 


and 










|H^(i27Tf)l 
^ w 


w. 


(45) 


(46) 


where 


H (s) 
u 


.004842 (1 + ,7425s) (1 + .2270s) 

(1 + 1.2569s) (1 + .134115) (1 +.2653s) 


(47) 


H (s) 
w 


.006331(1 + .7314s)(l+ ,11542s) 

(1 + 1.2410s) (1 + .06803s) (1 + .19733s)' 


(48) 


These linear-filter approximations to the power spectra of the 
rolling moments are also shown in figures 14 and 15, and it can 
be seen that they compare very closely with the more exact 
curves of Eggleston and Diederich. 


FLIGHT SIMULATOR EXPERIMENT 


iryr.^ dcscribes the flight simulator experiment, 

aircraft simulated, the flight Simula- 
pilot task/ the pilot's experience/ and the 

The purpose of the experiment was 
t determine whether or not pilots are sensitive to the differ- 
patchiness between Gaussian and non-Gaussian turbulence 
or the case of a STOL aircraft in a landing approach configuration. 


Vehicle Simulated 

It was decided to simulate a STOL aircraft on approach 
because the low frequency effects of turbulence become more 
pronounced as airspeed decreases. A STOL approach through 
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turbulence is a very demanding task for a pilot, and any differ- 
ences between turbulence models should be the most noticeable 
for this configuration. 

The aircraft simulated was the de Havilland DHC-6 Twin 
Otter, Series 300 (fig. 16) , chosen because its linear char- 
acteristics are fairly well documented, and because some flight 
time in an actual Twin Otter was available in order to aid the 
pilots in their assessment of the simulator fidelity. 

The configuration chosen was the following: STOL landing 

approach at a weight of 4989.5 kg., 36 m/sec., 40 of flaps, 
and center of gravity at 20% m.a.c. Data were available for 
a power-off approach configuration, but since a constant— heig 
task was desired (see page 45) and there were no data for this 
case, it was decided to use all the stability derivatives for 
the power-off case and simply add sufficient throttle to allow 
the airplane to fly level. No attempt was made to adjust the 
stability derivatives for the effects of power, other than by 
adding drag, lift, and pitching moment due to throttle changes 
from the reference flight condition. The aircraft is described 
in Table 3., and the equations used are discussed in Appendix B. 


Simulator 

The NASA-Ames 6 degree-of-f reedom, moving-base simulator 
was used (fig. 17). Because of mechanical stops, the cab 
travel is quite restricted (± 2.7 m (± 9 ft.) ) in the trans- 
lational degrees of freedom, but quite adequate in angular 
rotations (± 45° in roll, pitch and yaw) . Because of the 
tricted translational motions, it was necessary to wash out low- 
frequency accelerations with a high-pass filter in order to pre- 
vent them from integrating into motions which would cause the 
simulator to hit the travel limits. In addition, the simulator 
response was sluggish for high-frequency input signals, so 
compensation (in the form of a lead circuit) was added in order 
to improve the response in this frequency range. It was also 
found that large-amplitude input accelerations in the transla- 
tional modes caused the simulator drive system to shut off, so 
commanded accelerations from the equations of motion were put 
through a limiting circuit. Details of the washout, compensa- 
tion, and acceleration limits used are given in Appendix C. 

The instrumentation and controls were somewhat different 
between the simulator and the actual Twin Otter; however the 
pilots did not seem to feel these differences would have too 
much effect on their evaluation of the realism of the turbu- 
lence models. The simulator had a stick instead of control 
yoke, and a single throttle (for both engines) on a console on 
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TABLE 3. - AIRPLANE DATA 

Aircraft; de Havilland DHC-6 Twin Otter (Series 300) 
de Havilland data used where available. Missing derivatives 
were supplied by the authors. Some derivatives were altered 
for the purposes of this study. 

Airplane Parameters 


u = 36.03 m/sec (118.2 ft/sec) = 0 

U IT XLi 


ra = 4989.5 kg (341.9 slugs) 
S = 39.019 m^ (420 ft^) 
b = 19.812 m (65 ft) 

= 7.62 m (25 ft) 
c = 1.9812 m (6.5 ft) 

= 2.4384 m (8 ft) 


Ixx ^ 21,621 kgm^ (15,947 slug ft^) 
lyY = 31,824 kgm^ (23,472 slug ft^) 
1^2 = 48,857 kgm^ (36,035 slug ft^) 

Ijj 2 = 1.482 kgm^ (1,093 slug ft^) 
P__^ = 4.847X10^w (650 ff)per engine 

lUaX 

6t = .725 
0 

p = 1.2256 kg/m^ (.002378 slugs/ft^) 


Stability Derivatives 



0 


0.320 

^ a 

— 

- 0 . 090/rad 

a 

= 

0 . 9832/rad 

r 

= 

0.336 

0 

“ 

0.320 



- 0.777 

CT 

— 

- 0.75 

6a 


. 150/rad 

6 

— 

- 0.775/rad 

^5r 

= 

- 0.045/rad 

r 

— 

0.513 

CT 


- 0.1698 

P 


- 0.131 

C 

m 

a 

~ 

- 2.026/rad 

6a 

— 

0.0108/rad 

C 

m* 

a 


- 8.663/rad 

6r 

— 

- 0.391/rad 

^m 

q 

— 

- 28.76 

t3 

— 

- 0.39/rad 

c 

m » 
6e 

= 

2 . 068/rad 

0 

— 

1.576 

C 

m. 

tu 

— 

1.24 


42 



TABLE 3.- AIRPLANE DATA - Continued 


C = - 6.109/rad 

z 

a 

C = - 2.150/rad 

z • 

a 

C = - 7.144 

C = 0.5236/rad 

Z r 

6e 

C,. =0.55 

^CT 

C = - 0 . 632/rad 

^ta 

C = - 2.09/rad 

^ta 

C = - 2 .43/rad 


V- 

ta 

= - 8.04/rad 

C 

= 0.147/rad 

c 

nr 

= - 0.219 

C 

= - 0.0219/rad 

C 

= 0.1565/rad 

C 

^t3 

= 0.14/rad 

c 

= - 0.0935 


Resulting Characteristics of Airplane Dynamics 

a. ) Short-Period Mode 

natural frequency = 2.322 rad/sec 
damping ratio = .6538 

frequency of oscillation = 1.757 rad/sec 
period = 3.576 sec 

time to half-amplitude = .4566 sec 

b. ) Phugoid Mode 

natural frequency = .3265 rad/sec 
damping ratio = .2186 

frequency of oscillation = .3186 rad/sec 
period = 19.72 sec 

time to half-amplitude = 9.711 sec 

c. ) Roll Convergence 

time to half-amplitude = .1167 sec 
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TABLE 3.- AIRPLANE DATA - Concluded 


d. ) Spiral Divergence 

time to double-amplitude = 11.92 sec 

e. ) Dutch Roll 

natural frequency = 1.519 rad/sec 
damping ratio = .3431 

frequency of oscillation = 1.427 rad/sec 
period = 4.403 sec 

time to half-amplitude = 1.330 sec 
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the left side as opposed to two throttles overhead on the pilot's 
right in the airplane. The Twin Otter had the normal ILS cross- 
pointer, needle and ball, directional gyro, and artificial hori- 
zon. The simulator however had all this information on one 
instrument, the attitude director indicator. 


Description of the Experiment 

Turbulence Cases .- The tests were conducted for four cases 
composed of one set of measured turbulence data and three differ- 
ent statistical turbulence models, as follows. 

(1) Real - using time histories of real atmospheric turbu- 
lence recorded during the LO-LOCAT (low level critical air 
turbulence) program. 

(2) Gaussian model - a model having a Gaussian probability 
distribution. 

(3) Matched non-Gaussian model - a non-Gaussian model 
characterized by a standard deviation ratio of unity (see 
equation (39) and related discussion) . The probability 
characteristics of this model were chosen so as to be 

as much like those of real turbulence as possible. These 
characteristics are shown in figures 10 and 11 by the curves 
marked R = 1 . 

(4) More patchy non-Gaussian model - a non-Gaussian model 
characterized by a standard deviation ratio of two. The 
probability characteristics of this model were chosen so 
as to produce more large gusts and more patchiness than 
are found in real turbulence. The curves marked R = 2 
in figures 10 and 11 show the probability characteristics 
of this model. 

The real data and the three models were scaled to the same 
rms intensities, had the same power spectra and cross-spectra, 
and had the effects of spatial distribution included as outlined 
on page 35. The only difference between the models was their 
different probability densities, resulting in differing amounts 
of "patchiness". 

Pilot Task. - It was decided to have the pilots fly a cons- 
stant - altitude tracking task in order not to introduce too 
many variables (such as changing the turbulence characteristics 
with height) that might distract the pilots from their primary 
goal of trying to distinguish differences between turbulence 
models. This hypothetical task consisted of tracking a constant- 
width ILS localizer and a constant-altitude glideslope parallel 
to the ground at an altitude of 76.3 m (250 ft.). The sensitivity 
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of the ILS needles was set so that full-scale deflections corre- 
sponded to 15.2 m (50 ft.) off the glideslope and 76.3 m (250 ft.) 
from the localizer. These full scale deviations would occur at 
a range of 1750 m (5730 ft.) from the runway for a 2 1/2° glide- 
slope.) The pilot was told that his ILS tracking errors both 
in the vertical and horizontal were being measured, and that 
his primary task was to minimize these errors using rudder, 
aileron, elevator, throttle, and the flying techniques he would 
normally use on an approach. 

Test Procedures .- Four to six runs, each of five minutes 
duration, were made in one pilot session, and in the course of 
the experiment each pilot flew each of the turbulence models 
from three to five times. The models were given to the pilot in 
random order so that no consecutive runs had the same model. 
Steady crosswinds of either 2.5 m/sec (5 knots), - 3.8 m/sec 
(- 7 1/2 knots) , or 0 were fed in, and the pilot was told what 
the crosswind was. (He would receive this information from the 
tower in an actual approach.) After each run he was asked for 
his comments on the turbulence by means of a questionnaire (see 
Appendix D) . This question sheet asked the pilot to estimate 
the turbulence intensity and realism, his workload, his task 
performance, and a Cooper rating (ref. 23) for the airplane/ 
turbulence/ task combination. Additional questions tried to find 
out why the pilot did or did not find the turbulence realistic. 
These included comments on the relative amplitudes of distur- 
bances in roll, pitch, etc., whether the turbulence was too con- 
tinuous or too patchy, and whether the motions contained annoying 
amounts of high or low frequencies. 

During the flight a total of twenty-five flight parameters, 
such as angular rates and ILS errors, were recorded on magnetic 
tape for later analysis. 

Pilots' Experience .- Two pilots were used in the program, 
and their previous experience is summarized on the next page. 
Pilot A's 3 hours consisted of two flights in the three-week 
period prior to the experiment. These flights consisted of 
some level flight in smooth air at 36 m/sec and 40° of flaps 
to familiarize the pilot with the Twin Otter's stability and 
control characteristics, and in addition, some IFR landing 
approaches at 36 m/sec in moderate turbulence. Pilot B's 1.5 
hours were flown one year earlier. 
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PILOTS' EXPERIENCE 


PILOT 


TOTAL 
NO. OF 
HOURS 

HOURS 

ACTUAL 

INSTRUMENT 

FLYING 

HOURS IN 
SIMULATORS 

MAIN 
TYPES OF 
AIRCRAFT 
FLOWN 

HOURS 
IN TWIN 
OTTER 

Pilot 

A 

2800 

45 

140 

light twin, 
light single- 
engine, rotary 
wing 

3.0 

Pilot 

B 

4000 

100 

500 

Convair 340, 
light single- 
engine, Lear 
jet 

1.5 


RESULTS 


Results of Turbulence Simulation Development 

The development of a non-Gaussian turbulence simulation has 
been described under "Description of Atmospheric Turbulence" and 
in Appendix A of this report. It has been shown that, in theory, 
this new method can closely approximate the patchy nature, non- 
Gaussian probability distribution, and frequency content of 
atmospheric turbulence. Furthermore these properties can be 
varied during operation of the simulation in order to represent 
changing flight conditions such as vehicle altitude or airspeed. 

Digital computer programs modeling the system block dia- 
grams presented in Appendix A have been written and artificial 
turbulence time histories have been produced. A statistical 
analysis of these time histories has been carried out and the 
experimentally obtained results of the analysis of these time 
histories will now be compared with those expected from the 
theoretical considerations of Appendix A. The statistical pro- 
perties of interest in this report are; power spectral density, 
cross spectral density, probability distribution, and patchiness. 
Each of these quantities will now be considered in turn. Since 
similar results were obtained for all three statistical models 
considered, only the matched non-Gaussian model data will be 
presented . 

Power Spectral Densities .- Figures 18, 19, and 20 compare 
the power spectral densities of the three gust components at 
the vehicle center of gravity with the modified Dryden spectra 
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of figure 3. The fit of the theoretical curves is very good, 
especially in view of the fact that only one 5-minute data sample 
was used to produce the experimental curves. Because of the 
short time history length, spectral estimates below .1 hertz 
are quite uncertain. 

Cross Spectral Densities. - Figure 21 compares the desired 
cross spectral density of the u and w gust components with the 
results estimated from a 5-minute sample. The cross spectrum 
exhibits correct behavior for frequencies above .6 hertz. The 
large discrepancies below this frequency are attributed to 
uncertainty caused by the short time history length. 

Probability Distribution . - Figures 22 and 23 show typical 
smoothed cumulative probability and probability density distri - 
butions obtained from the simulated turbulence. Both figures 
show good agreement with theory. 

Table 4 compares the estimated fourth and sixth normalized 
central moments obtained from the analysis of the matched non- 
Gaussian model time history with the theoretically predicted 
values. 

TABLE 4.- COMPARISON OF THEORETICAL VALUES OF M. AND WITH 

4 6 

THOSE FROM ANALYSIS OF THE MATCHED NON GAUSSIAN MODEL TIME 

HISTORY 



^4 

^6 

From Analysis of Time History 

4.51 

42.73 

Theory (Eq. A-33) 

4.50 

52.50 


Agreement is excellent for M^ and fair for Mg . The error in 

Mg is attributed to the normal difficulty in estimating the 

higher moments of a random function, and to possible non- 
Gaussian behavior of the random number generator; it is not 
thought to indicate a defect in the theory. 

Patchy Character .- There is no theoretical model for the 
patchy nature of turbulence which is to be reproduced by the 
turbulence simulation, therefore no comparison with theory can 
be presented here. Instead, the reader's attention is directed 
to figures 4 and 13, which demonstrate the ability of the non- 
Gaussian turbulence simulation to model the patchy character 
of atmospheric turbulence. 
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Results of Flight Simulator Experiment 

Simu lator Fidelity ." The resemblance of the simulator fly- 
ing characteristics to those of the Twin Otter was considered 
to be quite close. One pilot thought the roll sensitivity to 
aileron might be a bit low, however had already been 

6a 

increased by 50% from the de Havilland data, and the other 
pilot was happy with it. 


The main criticism of the simulation was the presence of 
a lot of high-frequency noise, which showed up as a shudder or 
vibration when simulating smooth air flight. This made the 
response tcp turbulence seem quite sharp-edged, as opposed to 
the wallowing motion felt in the Twin Otter. It is felt that 
two sources are to blame for this. First and most important, 
the fact that these vibrations occurred for smooth air condi- 
suggests that the simulator— drive system and simulator 
structure, coupled with the compensation network used to improve 
the high-frequency response (see Appendix C) were to blame. 

A second factor is probably the neglect of unsteady-aerodynamic 
effects of the gusts acting at the c.g. and the tail. A sharp 
gust at the c.g. will immediately become a similar vertical 
acceleration through the C stability derivative. In real 

* ^ ct 

the unsteady effects require a finite time for the lift 

to build up on an aerodynamic surface subjected to an instan- 
taneous change in angle of attack. The pilots felt that this 
unwanted noise, whatever the source, was sufficient to mask 
the small differences that the pilots were trying to detect 
between the various turbulence models. 

Another complaint which the pilots felt detracted from 
the realism was the absence of any drifting in the mean wind. 
Since the crosswind was set to a constant value for the whole 
flight, the pilot would just set his heading to correct for 
this and then wouldn't have to worry about it any more. On a 
approach, a pilot flies through a crosswind of varying 
intensity as he changes altitude, consequently he has to be 
continually hunting for a bias heading and a bias power-setting 
otder to stay on the localizer and glideslope. The turbu- 
lence model used has the capability of changing the turbulence 
characteristics with altitude, but it was felt that since the 
task was a constant-altitude one, the gust parameters should 
be fixed at those of the altitude flown. Obviously this 
detracted from the realism. 


Results of Pilot Questionnaire .- Answers from the qqestion 
sheets are shown in figures 24 to 30, and 32 to 34, but for 
brevity not all the answers are discussed here or shown in these 
figures. The only plots shown here are those of special interest 
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or those for which there was a difference between answers for 
the various turbulence models. Results for the remaining ques- 
tions are shown in Appendix E . 

These plots of the pilots' answers to the questionnaire 
are to be interpreted as follows. Each individual answer is 
shown by a small "x" or square. An average of the points for 
each pilot and model is shown by the large "X" or square. Using 
the Student's t distribution (ref. 24), and (for the case of 
turbulence intensity, figure 24) assigning the value 1.0 to the 
answer "light", 2.0 to "moderate", and 3.0 to "severe", one can 
calculate the 95% confidence limits for where the true mean must 
lie for a given pilot and mode. These limits are shown by the 
vertical bars in the figures. In other words, the data points 
have a mean denoted by the large "X" or square. Based on these 
data points, there is 95% confidence that the true mean for a 
given pilot and model lies somewhere in the region denoted by 
the vertical bar. 

These confidence limits for where the true mean lies give 
a way to tell whether or not a pilot can readily tell any differ- 
ence between two turbulence models. One can be 95% confident 
that the pilot can distinguish between two models if the follow- 
ing two conditions are both met: 

(1) The mean for the first model must lie outside the 

95% confidence limit for the mean of the second model; and 

(2) The mean of the second model must lie outside the 

95% confidence limit for the mean of the first model. 

An example where the pilot can differentiate (at the 95% confi- 
dence level) between models is pilot B's ratings of Gaussian 
and Real turbulence in figure 24. An example of the pilot being 
unable to distinguish between two models is shown in figure 28 
for pilot A's ratings of Gaussian and Real turbulence. In this 
case, the mean for Real is outside the confidence limit for 
Gaussian; however the mean for Gaussian is within the confidence 
limits for the Real mean. 

Turbulence Intensity .- Inspection of the pilot estimates 
of turbulence intensity shown in figure 24 shows that both 
pilots were able to distinguish between the Real turbulence 
and the Gaussian model - the Real one appearing to be the more 
intense. This is interesting, since the models were all scaled 
to have the same rms values. Apparently, because of its more 
continuous nature, the Gaussian turbulence seemed less severe. 

The pilots were unable to distinguish between any of the other 
models on the basis of turbulence intensity. 
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Realism of Turbulence .- Estimates for the turbulence realism 
are plotted in figure 25. Pilot A found no differences, rating 
all the models "fair" on the average. Pilot B was able to dis- 
tinguish between the Real and Gaussian models and between the 
Real and matched non-Gaussian ones. His mean for the Real turbu- 
lence was the highest, being slightly better than "good". 

Correctness of Relative Amplitudes of Disturbances .- Plots 
for the opinion of correctness of roll and yaw amplitudes are 
shown in figures 26 and 27 respectively. Pilot B generally 
seemed to think the amplitudes were about right, whereas pilot 
A (who had the recent Twin Otter time) seemed to feel there 
was a bit too much yaw and not quite enough roll in all the 
models. The pilots were happy with the amplitudes of the dis- 
turbances in the other axes, except that pilot B thought there 
was not enough side-force. 

Patchy Characteristics .- The patchiness of the turbulence 
is the one characteristic which is different in the four cases. 
The Gaussian model turbulence is quite monotonous or continuous 
in its amplitude, whereas in nature the turbulence tends to have 
periods of relative quiet followed by relatively noisy bursts. 

It was hoped that the pilots could distinguish between "patchy" 
models and a relatively continuous one. The pilots' answers to 
the question on patchiness are plotted in figure 28. It was 
hoped the pilots would rate the Gaussian model too continuous, 
the Real turbulence and the matched non-Gaussian model about 
right, and the more patchy non-Gaussian perhaps too patchy. As 
the figure shows, this was not quite the case, although the 
pilots were able to distinguish between some of the models. The 
Gaussian one was found to be "a little too continuous" , and the 
other models all bracket being "about right", although their 
mean values all tend to be slightly towards being continuous. 

Statistically, pilot A was able to distinguish that the 
matched, non-Gaussian and the more patchy non-Gaussian models 
are patchier than the Gaussian. He was unable to distinguish 
between the Real and Gaussian, since the Real has such a wide 
spread on the confidence limits for the mean. Pilot B could 
tell that the Real and more patchy non-Gaussian models were 
patchier than the Gaussian. He was unable to distinguish bet- 
ween any of the other cases. 

Frequency Content of Turbulence .- Pilot estimates of low 
and high-frequency content are shown in figures 29 and 30. Both 
A and B rated the low frequency content as being somewhere bet- 
ween "about right" and "not enough". Although pilot B rated 
nine of the flights as being about right, he continually 
remarked that the models lacked what he called "low- low fre- 
quency content" , which would cause him to have to hunt for the 
correct directional bias required to stay on the localizer. 
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(See the discussion on page 49.) For the high frequency criteria, 
pilot A consistently complained of too much high frequency noise 
from the simulator drive system. Pilot B, who did not have any 
recent Twin Otter time, rated the high frequencies as about 
right. (Pilot B has had considerably more simulator time, and 
might be more prone to filter out simulator shortcomings.) 

However a third pilot, who took part in the check-out and who 
had several hours of recent Twin-Otter time, also complained 
of too much high-frequency noise. 

Cooper Rating .- The pilots were asked to rate the airplane/ 
turbulence/ task combination using the Cooper-Harper Pilot Rating 
Scale shown in figure 31. In figure 32, pilot A's mean ratings 
were in the 4 1/2 to 5 range, while B rated the combination about 
a full point lower. Pilot A's means show no statistical differ- 
ence between models, whereas B rated the Real case higher than 
the Gaussian and the matched non-Gaussian. This corresponds to 
the higher rating for turbulence intensity (figure 24) found 
for the Real turbulence. 

Root Mean Square Values of Airplane Motions .- Root-mean- 
square (rms) values were computed for some of the airplane 
motions to see if these showed any differences between models. 

Rms values of the aircraft attitude angles 0 , (p , and 
are plotted in Appendix E (figures E-6, E-7, and E-8) and show 
no consistent differences between the turbulence models. How- 
ever the rms values of , the height above the glideslope, 

shows an interesting result in figure 33. Both pilots show a 
statistically higher rms value for the Real turbulence than 
for the Gaussian. This fits in with the estimates of higher 
intensity for the Real than for Gaussian turbulence, and is 
significant because the models had the same gust rms values. 

In comparison, the rms values of displacement from the localizer 
(fig. 34) show no significant variation between gust cases. 


CONCLUSIONS 


This report has described a new method for producing arti- 
ficial turbulence time histories which is more flexible in its 
ability to match desired statistical properties than are pre- 
sently used methods. The use of this method should give 
improved realism and accuracy when in piloted simulator studies, 
for aircraft response studies which investigate stability, 
stability augmentation, guidance, and for studies of aircraft 
structures as affected by atmospheric turbulence. 
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The new method was used to produce time histories of turbu- 
lence, which were analyzed in order to obtain certain of their 
important statistical properties. These properties were com- 
pared with the properties of actual turbulence. 

The conclusions from the analytical studies of the time 
histories developed by this new method as compared to those 
developed by other simulation methods are outlined below. 

1. ) The new method is capable of producing a class of non- 
Gaussian probability density distributions, which has not 
been possible with previous turbulence simulations. Using 
this method one can produce probability distributions 
which closely match the various probability distributions 
which are obtained from the analysis of real turbulence 
data. 

2. ) The results of the analyses show that the patchy 
nature of real turbulence is matched by this method. 

The use of this method in piloted simulators should 
remove the objection that the use of the Gaussian pro- 
bability results in turbulence which is too monotonous. 

3. ) The frequency content of the real turbulence can be 
matched. Although other methods can do this as well, 
this important property of turbulence simulation is 
preserved in this method. 

4. ) Although not demonstrated experimentally this new 
method has the capability common to other simulation 
techniques of simulating the time-varying characteristics. 

5. ) Since the probability density distribution can be 
closely matched to that of real turbulence, this method 
should be capable of predicting the occurrence of large 
gust velocities more accurately than present methods. 

This should prove especially valuable in applications 
in the structures area. 

6. ) While the mechanization of the proposed method is 
not as simple as that of the Gaussian model, it is still 
easily within the capabilities of modern digital or 
analog computers. 

The time histories derived from the non-Gaussian statis- 
tical models along with those from real turbulence and from 
a Gaussian model, were used in a flight simulator experiment 
in order to compare the realism of these various methods for 
simulating turbulence. The principal conclusions resulting 
from the flight simulator experiment are as follows. 
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1. ) The pilots tended to rate the Gaussian model as being 
"a little to continuous," whereas the two non-Gaussian 
models and the real time history were rated "about right" 
with regard to patchiness. 

2. ) The pilots comparison of intensity indicated that, 
even though all models were scaled to the same intensi- 
ties, the non-Gaussian models seemed to be more severe 
than the Gaussian model. It is concluded that this is 
due to the fact that there are more large gusts when 
using the non-Gaussian models. The rms of the glide 
path errors reinforced this view, as these errors were 
definitely greater for the real turbulence and the non- 
Gaussian models than for the Gaussian model. 

3. ) Although the pilots expressed the impressions discussed 
above, other impressions and statistical analyses of pilots' 
performance were not conclusive in differentiating among 
the various turbulence models. One possible reason for 
this might be that, even though the non-Gaussian model 
extended the range of statistical properties matched, 
pilots may be sensitive to certain properties which are 

not well matched by this model. A second reason might be 
that the characteristics of the simulator itself, mainly 
its high-frequency noise characteristics, tended to mask 
some of the differences in the turbulence models. Another 
possible reason might be that the pilots had limited 
flying time in turbulence in the actual aircraft being 
simulated . 


Department of Aeronautics and Astronautics, 
University of Washington, 

Seattle, Washington, March 30, 1973. 
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FIGURE 3. NORMALIZED GUST POWER 
SPECTRAL DENSITIES ASSUMED IN 
THIS REPORT 
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FIGURE 5- LINEARLY-FILTERED GAUSSIAN WHITE-NOISE 
SIMULATION 
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FIGURE 6. PATCHY CHARACTERISTIC OF ACTUAL 
ATMOSPHERIC TURBULENCE COMPARED WITH THAT 
PRODUCED BY THE GAUSSIAN TURBULENCE MODEL 
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FIGURE 7. COMPARISON OF MODIFIED 
BESSEL FUNCTION DISTRIBUTION WITH 
GAUSSIAN AND MEASURED DISTRIBUTIONS 
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FIGURE 8. PATCHY CHARACTERISTICS OF MODIFIED BESSEL 
FUNCTION MODEL COMPARED WITH THOSE PRODUCED BY THE 
GAUSSIAN MODEL AND ACTUAL ATMOSPHERIC TURBULENCE 


NaissnbQ 


cn 

CEljJ 
lUh- 
Z _J 



UJ 

CO 

•-I CO 
OLU 
Z U 
I QC 
UJZ5 
1—0 




65 


FIGURE 9- PROPOSED NON-GflUSSION TURBULENCE SIMULnilON 
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FIGURE 10. COMPARISON OF CUMULATIVE 
PROBABILITY DISTRIBUTIONS FOR VARIOUS 
STANDARD DEVIATION RATIOS 





FIGURE 11. COMPARISON OF PROBABILITY 
DENSITIES FOR VARIOUS VALUES OF THE 
STANDARD DEVIATION RATIO 
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FIC5URE 12. DEPENDENCE OF THE FOURTH 
AND SIXTH NORMAL I ZEO CENTRAL MOMENTS 
UPON THE STANDARD DEVIATION RATIO 
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FIGURE 13. DERIVATIVE OF SIMULATED TURBULENCE 
SHOWING DEPENDENCE OF PATCHY CHARACTERISTIC 
UPON THE STANDARD DEVIATION RATIO R 
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FIGURE 14. POWER SPECTRUM OF ROLLING 
MOMENT DUE TO DISTRIBUTION OF HERD-ON 
GUSTS ACROSS SPAN 
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FIGURE 15. POWER SPECTRUM OF ROLLING 
MOMENT DUE TO DISTRIBUTION OF VERTICAL 

GUSTS ACROSS SPAN 
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FIGURE 17. NflSfl-fltlES 6-DEGREE FLIGHT SIMULPTOR 
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FIGURE 19- POWER SPECTRAL DENSITY 
OF EXPERIMENTALLY GENERATED LATERAL 

GUST COMPONENT 
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FIGURE 21. CROSS SPECTRPL DENSITY OF EXPERIMENTALLY 
GENERATED GUSTS (MATCHED NON-GAUSSIAN MODEL) 




CUMULfiTlVE PROBFIBILITY.P(^) 



NORMALIZED GUST VELOCITY.^ 
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GENERATED TURBULENCE 
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FIGURE 23. PROBABILITY DENSITY 
DISTRIBUTION OF EXPERIMENTALLY 
GENERATED TURBULENCE 
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FIGURE 31. COOPER-HARPER PILOT RATING SCALE 
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APPENDIX A 


DETAILED DESCRIPTION OF TURBULENCE SIMULATION 


The purpose of this appendix is to describe in detail the 
turbulence simulation technique proposed in the main body of 
the preceding report. This discussion will be divided into 
six sections. 

A brief review of the statistical properties which are 
to be simulated. 

Discussion of turbulence simulation techniques now avail- 
able and arguments leading to the technique proposed in 
this report. 

A mathematical analysis of the system which generates the 
uncorrelated components of the turbulence simulation. 

Analysis of correlation techniques by which the indivi- 
dual components of the turbulence simulation can be 
related. 

Description of a coordinate rotation scheme by which the 
specific cross spectral densities needed for the flight 
simulator experiment of this report were produced. 

Construction of a complete block diagram of the proposed 
turbulence simulation, and specialization of this general 
system to the three artificial turbulence models used in 
the flight simulator experiment of this report. 

This appendix is intended to function as a self contained 
description of the proposed non-Gaussian turbulence simulation. 
It therefore contains symbols, illustrations, and references 
distinct from those of the preceding report. 
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SYMBOLS 


a 

b 


c 

C. 


'K 


ug 


wg 


uw 

d 

erf 


E 

f 


F 

9 - 


h. 

J 

H . 
3 


j 

K . 
3 


K 


o 


random time history 
random time history 
random time history 
Gaussian correlation parameter 

modified Bessel Function correlation parameter 

non-dimensional rolling moment induced by longitudinal 
gust component distributed along wing span 

non-dimensional rolling moment induced by vertical 
gust component distributed along wing span 

cross correlation of u and w gust components 

random time history 

error function 

expected value operator 

frequency (hertz) 

Fourier transform operator 

Gj (s) 

(s) 

source n j 

modified Bessel function of second type, order zero 


impulse response of linear filter 
linear filter with transfer function 
impulse response of linear filter 
linear filter with transfer function 

integer 

power spectral density of white noise 

2 

in units of (output variable) /hertz 
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scale length of the gust component ( m ) 

normalized central moment of a random process 
defined by the equation 

1 ,x(t)-(mean value of x) , ^ 

M. = lim ^ ^ i 

J T^-oo -T X 


integer 

Gaussian white noise process 
probability density distribution 
cumulative probability distribution 

ratio of modified Bessel process standard deviation 
to Gaussian process standard deviation (see eq. A— 30) 

Laplace variable 

time (sec) 

gust component in direction of mean wind relative 
to surface ( m /sec) / also denotes simulated turbu 
lence time history ‘ 

longitudinal gust component in vehicle coordinates 
( m /sec) 

mean true airspeed of vehicle ( m /sec) 

lateral gust component in earth-fixed coordinates 
( m /sec) forming right hand system with u and w 
( m /sec) 

lateral gust component in vehicle coordinates ( m /sec 

vertical gust component in earth— fixed coordinates 
positive downward ( m/sec) 

vertical gust component in vehicle coordinates 
( m /sec) 


random time history 



y random time history 

z random time history 

Y variable of integration 

6 Dirac delta function 

n Gaussian white noise 

A variable of integration 

5 variable of integration 

standard deviation of the random process x(t) 

(units of x) 

T correlation time increment (sec) 

(|) power spectral density of the random function x(t) 

^ 2 

(units of X /hertz) , <t> is an even function of 
frequency, related to a by the equation 

cx> 

o ^ = f 4) (f) df 

X ^x 

— oo 

4> cross spectral density of random processes X(t) and 

^ Y(t) (units of XY/hertz), <() is the Fourier transform 

of the cross correlation of X and Y 

^ heading angle of vehicle relative to the mean wind 

in earth-fixed coordinates 

( ) binomial coefficient 

Operators: 

* convolution (see eq. A-45) , the free variable of the 

convolution is denoted by a subscript 

_ (underbar) complex conjugate 


90 



Statistical Properties to be Simulated 


The turbulence simulation described here consists of a 
linear filter array driven by Gaussian white noise generators. 
The output of this system consists of five simultaneously 
generated random time histories which are to be inputs to a 
flight simulator. These five time histories can be physically 
interpreted as follows. 

1. ) The longitudinal gust component (in earth-fixed 
coordinates, parallel to the mean wind) occurring at 
the vehicle center of gravity. 

2. ) The lateral gust component (in earth-fixed coordinates) 
occurring at the vehicle center of gravity. 

3. ) The vertical gust component (in earth— fixed coordi- 
nates) occurring at the vehicle center of gravity. 

4. ) The non-dimensional rolling moment of the vehicle 
(in body axes) due to the distribution of longitudinal 
gusts along the wing. 

5. ) The non-dimensional rolling moment of the vehicle 
(in body axes) due to the distribution of vertical 
gusts along the wing . 

In addition to these five time histories, the vehicle yawing 
moments and gusts occurring at the tail surfaces are also 
reguired for a realistic flight simulation. Yawing moments 
due to the distribution of gusts along the wing are assumed 
to be directly proportional to the rolling moments (inputs 
4 and 5 above) and therefore need not be explicitly generated 
by the turbulence simulation. Furthermore, it is assumed 
that some means is available to delay the gust components 
(inputs 1, 2, and 3 above) in order to provide gusts at the 
vehicle tail surfaces. 

The time histories produced by the turbulence simulation 
are to model turbulence and its effects in a statistical sense. 
The statistical properties of interest in this report are the 
following. 

1. ) Power spectral density, the average frequency content. 

2. ) Cross spectral density, the average phase relationship 
between the frequency components of two time histories. 
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3. ) Probability distribution, a description of the values 
assumed by a random time history and the probability of 
their occurrence. 

4. ) Patchiness, the tendency of turbulence to occur as 
motion^ activity separated by regions of less intense 

Each of these four properties will now be discussed in detail. 

Power Spectral De nsity of Artificial Tnr hni^n^o _ The 
spectral densities of the gusts at the vehicle center 
gravity are a modified form of those suggested by H L 
Dryden (reference A-1) . y . j-.. 
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(A-3) 


that the form of 

^ pectrum has been changed to a form similar to 

longitudinal component, and independent 
scale lengths have been assumed for each component. It is 
important to note that Dryden 's equations were intended to 

only stationary, isotropic, homogeneous, and ergodic 
turbulence. The equations represented above are, in general, 
intended to represent non-stationary , non-iso tropic, and non- 
homogeneous turbulence. In order to allow this generality 
the scale lengths, L , and standard deviations, a , of equa- 
tions (A-1) through (A-3) must be functions of altitude, 
atoospheric stability, vehicle heading relative to the mean 

velocity, and any other variables which may 
ect the power spectra. Since the turbulence considered in 
this report is assumed to be stationary, homogeneous and 

^ constants in the particular 

turbulence models described here. ^ 
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The power spectral densities of the rolling moments induced 
by the distribution of gusts along the wing are dependent upon 
the geometry of the vehicle as well as the characteristics of 
turbulence. Different spectra will be required for each vehicle 
to be simulated. (A technique for deriving these spectra is 
presented in reference A-2.) The particular spectra suitable 
for simulating the DHC-6 "Twin Otter" aircraft flying at an 
altitude of 76.2 m (250 ft.) above the surface in the landing 
configuration is 

a ^ (.004842) ^ [1 + .5513 (2Trf)^J [1 + . 0515 (2TTf ) ^] 

. = u ^ 2 

[1 + 1.5798 (2TTf) [1 + .0180 (2irf) [1 + .0704 (27Tf) ] 

ug 

a ^ (.006331) ^[1 + .5349(2TTf)^] [1+ . 0133 (2TTf ) ^ ] 

. ^ w 2~2 

[1 + 1.5401 (2TTf) [1 + .004628 (2-rrf) [1 + .03894 (2TTf) ] 

wg 


The dependence of these spectra upon the gust scale length and 
the vehicle airspeed has not been given as in equations (A-1) 
through (A-3) because these spectra were sufficient for the 
flight simulator experiment described in this report and further 
generalization was considered unnecessary. The simulation of 
a different vehicle or different operating conditions will 
require the derivation of new rolling moment spectra. 

Cross Spectral Densities.— The cross correlation of gust 
components is typically neglected in turbulence simulations 
because explicit forms for cross spectra are usually not known 
and because the generation of cross spectra greatly complicates 
the simulation process. Cross spectra were considered necessary 
for the flight simulator experiment described in the main body 
of this report because its object was to compare the realism of 
turbulence simulations which differed only in their probability 
distributions and patchy characteristics. The matching of cross 
spectra was necessary in order to eliminate the correlation of 
turbulence components as a variable of the experiment. 

The turbulence simulation described here produces a non- 
zero cross spectrum between the longitudinal and vertical gust 
components at the vehicle center of gravity. All other cross 
spectra are identically zero. The decision to correlate only 
the longitudinal and vertical gusts is based upon evidence pre- 
sented by Elderkin (reference A-3) which indicates that only 
these components are correlated. It should be noted that 
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Elderkin's data were obtained from tower measurements, 
ent results might be obtained from airborne measurements 
while flying across the mean wind. 


Differ- 

taken 


The cross spectral density chosen to relate the u and w 
gust components is based upon only a qualitative analysis of the 
data presented in reference A-3. The mathematical form is suq- 
gested by the method of cross correlation used. 


<u„> 


[1+3 (TTCj^f ) -i (2TTCj^f ) ] 
[l+(TTCj^f)^]^ 
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(A- 6) 


where R , , and are arbitrary parameters satisfying 

the inequalities. 


C„ > 




R 

> 0 


(A-7) 

2L 

w 

u 

0 


^0 

(A-8) 

L 

w 


L 

> 

u„ 

(A-9) 


The physical significance of these parameters in the simulation 
system will be described in a later section of this appendix. 


Probability Distribution .- In this report, 
characterizing the probability distributions 
will be used. 


three methods 
of turbulence 


1. ) Cumulative probability distribution. 

2. ) Probability density distribution. 

3. ) Normalized central moments. 


94 



Figures A-1 and A-2 present typical cumulative probability 
and probability density distributions of atmospheric turbulence 
from reference A— 4. In both cases the result expected of a 
Gaussian process is also given for comparison. Note that atmos- 
pheric turbulence has more large and small gusts than a Gaussian 
process . 

Another way to consider the probability distribution of a 
random process is through its moments. Table A-1 below compares 
the fourth and sixth normalized central moments of a typical 
turbulence time history with those expected from a Gaussian pro- 
cess. The higher values measured in real turbulence indicate an 
increased probability of large gusts. The real turbulence time 
history used in Table A-1 was a five minute sample from the 
LO-LOCAT data collection project (LO-LOCAT Phase 3, Test 1141). 

TABLE A-1.- COMPARISON OF THE FOURTH AND SIXTH NORMALIZED 
CENTRAL MOMENTS OF ATMOSPHERIC TURBULENCE WITH THOSE OF A 
GAUSSIAN PROCESS 



^4 

^6 

Atmospheric turbulence 

3.5 

21.7 

Gaussian process 

3.0 

15.0 


The data presented above indicate, at least qualitatively, 
the non-Gaussian behavior of atmospheric turbulence. In view 
of the limited data available, no explicit form for the prob- 
ability distribution will be suggested here. However, the turbu- 
lence simulation to be proposed in the next section of this 
appendix will be required to have a distribution with the 
general characteristics described above. 

So far nothing has been said about the probability distri- 
bution of vehicle rolling moments, or possible differences in 
the distributions of the three gust components. In view of the 
limited data describing the gust distributions and the complete 
absence of data describing the rolling moments, all components 
of the turbulence simulation will be assumed to have the same 
probability distribution. 

Patchy Characteristics .- There is no easily formulated 
statistical quantity representing the patchy nature of atmos- 
pheric turbulence. However, the existence of a patchy character 
can easily be observed in recordings of the gust time derivative. 
Figure A— 3 compares a typical gust derivative time history with 
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firoin two of ths ciirtif ici^l turbulsncs ^snsirstion 
methods to be considered presently. Note that the turbulence 
time history has a clearly discernable patchy nature. 

Not only should each component of the turbulence simulation 
have this characteristic, but examination of measured data 
shows that the patchiness is correlated for all three components, 
^his correlation of patches will be discussed further below. 


Turbulence Simulation Techniques 

. this section various methods of producinq time histories 

with the statistical properties described above will be considered. 

The Gaussia n Simulation .- By far the most widely used turbu- 
lence simulation is the linearly filtered Gaussian white noise 
method. A block diagram of this system is shown in figure A-4 . 

The technique is quite simple and consists of passing Gaussian 
white noise through a linear filter. The filter transfer func- 
is chosen so as to produce an output with the power spec- 
tral density of turbulence. Proper choice of the transfer func- 
tion allows virtually any power spectrum of practical interest 
to be produced. The spectra proposed in this report, equations 
through (A-5) , are easily generated. It is also possible 
P^o^uce cross spectra between two or more components of the 
simulated turbulence. Thus the linearly filtered Gaussian white 
noise turbulence simulation can be used to accurately model the 
frequency characteristics of turbulence. 

Unfortunately, this technique cannot reproduce the prob- 
abalistic nature of turbulence because the time histories which 
It produces are always Gaussian. As shown in figures A-1 and 
A-2 atmospheric turbulence is clearly non-Gaussian. Therefore 
the Gaussian model is unrealistic. The linearly filtered Gaussian 
white noise method also cannot model the patchy nature of turbu- 
lence, Figure A- 3 compares the patchy characteristics of a 
Gaussian process with those of real turbulence. Note the complete 
absence of patches in the Gaussian time history. The frequency 
content of both time histories is the same, the different char- 
acter of the Gaussian process is due entirely to the fact that 
it is Gaussian. 


The M odified Bessel Simulation .- The patchy nature of atmos- 
turbulence suggests that a simulation based upon the multi- 
plication of two Gaussian processes may provide a realistic repre- 
sentation. The first process would represent constant intensity 
turbulence and the second would represent its time varying stan- 
dard deviation. Figure A-5 presents a block diagram of such a 
system. 
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A turbulence simulation of this type has already been 
studied in reference A~5 where it is shown that this technique 
is capable of modeling the frequency content of turbulence just 
as is the Gaussian model described above. From the standpoint 
of probability distribution and patchy characteristics, however, 
the product of Gaussian processes appears to be too extreme. 
Figures A-6 and A-7 present the cumulative probability and 
pi^obability density produced by this simulation. (The probabil~ 
ity density function is proportional to the modified Bessel 
function of second type and order zero, hence the name "modified 
Bessel simulation.") Figure A-3 compares the patchy characteris- 
tics to those of the Gaussian simulation and real turbulence. 

Notice that in all three figures (A-3, A-6, and A-7) the 
Gaussian and modified Bessel processes appear to "bracket the 
behavior of real turbulence. The modified Bessel simulation 
produces more large gusts, more small gusts, and more severe 
patches than are found in real turbulence while the Gaussian 
simulation underestimates each of these. It therefore seems 
reasonable to propose a turbulence simulation which combines 
both the Gaussian and modified Bessel simulations. Figure A-8 
shows the block diagram of such a system. 

The statistical properties produced by the system of figure 
A-8 will be developed in the next section of this appendix. 
However, certain characteristics are immediately apparent and 
will be discussed here. 

The random process c(t) in figure A-8 is produced by the 
same system shown in figure A-5; it is therefore a modified 
Bessel process. Similarly, d(t) is produced by the same 
system shown in figure A-4 and is therefore a Gaussian process. 

If the standard deviation of the Gaussian process is made much 
less than that of the modified Bessel process, u(t) will 
become essentially modified Bessel in nature. On the other 
hand, if the standard deviation of c(t) is made small rela- 
tive to that of d(t) , the system output will become Gaussian. 
It is clear, therefore, that the proposed turbulence simulation 
is capable of producing either the Gaussian or the modified 
Bessel simulations as special cases. In general the system 
output will lie somewhere between these extremes. 


Statistical Analysis of Turbulence Simulation 

A statistical analysis of one component of the proposed 
simulation (figure A-8) will new be presented. Expressions 
for the power spectral density and probability distributions 
of the system output will be derived- In figure A-8 recall 
that n, , Ho f and n, are statistically independent Gaussian 
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white noise sources, and , G 2 , and G 3 are linear filters 

with transfer functions Gj^(s) , G 2 (s) , and G^ (s) respectively. 

The time histories a{t) and b(t) are multiplied to produce 

c(t) which IS then summed with d(t) to give the overall system 
output. ^ 

Power Spectral Density .- In order to derive an expression 
for It will first be necessary to consider the power 

spectra of a(t) , b(t) , c (t) , and d(t) . The well known 
^fi^tionship between the input and output power spectral densi- 
ties of a linear filter with transfer function G(s) is 


$y(f) = |G(i27Tf) |2 $^(f) 


where y(t) is the system output and x(t) is the system input. 
It x(t) IS a white noise process then is a constant, 

independent of f . In this case the shape of $ is deter- 
mined entirely by the filter transfer function. In figure A- 8 , 
a , b , and d are the responses of linear systems to white 
noise inputs. Denote the power spectral constants of the white 
noise sources , 112 , and 113 by , K 2 , and K 3 respec- 

tively. Then by application of equation (A-10) the power spec- 
tral densities of a , b , and d are: 



|G^(i27Tf) 

l\ 

(A-11) 

♦b'*> = 

|G2(i27Tf) 


(A-12) 

t^(£) = 

|G3(i2TTf) 


(A-13) 


Next consider the power spectral density of c(t) . It 
is proven in reference A-5 that if z(t) is the product of 
statistically independent random processes x(t) and y(t) 
then the power spectral density of z (t) is 


^^(f) = f $ (f-X)dX (A-14) 

— 00 ■* 


In the present case c(t) is the product of a(t) and b(t) 
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Furthermore a(t) and b(t) are statistically independent because 
they are derived from statistically independent noise sources. 
Therefore the power spectral density of c(t) is 


$ (f) = / $ (X) (f-X)dX (A-15) 

c a b 


An expression for r the power spectral density of the 

system output, can now be given. The process u(t) is the 
sum of c(t) and d(t) . These two random time histories must 
be statistically independent because they are generated by 
statistically independent white noise sources. The power spec- 
tral density of a sum of independent random processes is merely 
the sum of their respective power spectra. Therefore 




(A-16) 


Substitution for $ from eq. (A-15) and for ' '^d 

from eq. (A-11) through (A-13) gives 

$ (f) = / iG, (i2TTX) |2|G,[i2iT(f-M] l^dX + K lG3(i27Tf) 1^ 

(A-17) 

Equation (A-17) is a general expression for the power spectral 
density of the system output. Considerable simplification 
results if the filter transfer functions r ^2 ' ^3 

required to satisfy the condition 


K3lG3(i2TTf) 


K1K2 
,, 2 „ 2 
°a ^b 


/ lGj^(i2TrX) |^|G2(i2Ti(f-X) 1^ dX 

(A-18) 


In this case the power spectral shape becomes invariant with 
respect to the standard deviations ‘ 

expression for becomes 


t„(f) = 


K3lG3(i2TTf) 


(A-19) 
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Note that the shape of is now determined entirely by the 

choice of just as in the linearly filtered Gaussian white 

noise turbulence simulation. 


In the following it will always be assumed that the condi- 
equation (A-18) is satisfied. Hence, equation 
(A-19) gives the power spectral density of the system output, 


Proba bility Distribution .- When the linear filters of the 
system have been chosen so that condition (A-18) is satisfied 
the shape of the output power spectral density is not influenced 
by changing the standard deviations of c (t) or d(t) . However 
the probability distribution of the output is affected. Expres- 
P . 3P , and M will now be derived. 


u 


1 .) 


u ' "n 

Probability Density Distribution, P 


u 


Consider again figure A-8. The white noise time histories 
f 1 I 2 » and are, by definition, Gaussian. It is well 

known that when a Gaussian process is passed through a linear 
system it remains Gaussian. Therefore a{t) , b(t) , and d(t) 
are Gaussian and the probability density of each is 'the familiar 
normal distribution. 


P^(x) = 
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(- 1 (^) 1 
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(A-20) 

0 /Ttt 

a 

exp 

P^(x) = 
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exp 
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(A-21) 


P^(x) = 

1 

exp 

1 2 

[- — (JL) ] 

d 

(A-22) 


P<-(X) . 

The 

time 

history c(t) 

is the pro- 


duct of a(t) with b(t) . A derivation of the probability 
density distribution of a product is presented in reference 
A-5. The result is 





(A-23) 
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In deriving an expression for P^(x) it is convenient to 
first find IP^(x) , then obtain P^(x) by differentiation. 

By definition 


3P^(x) = PROBABILITY [u (t) < x] 


(A-24) 


In this case 


u ( t) = c (t) + d (t) 


Therefore expressed symbolically as 


oo d (t) — ^ 

IP (x) = / PROBABILITY [AND SIMULTANEOUSLY] dX (A-25) 

U -00 c(t) <x-X 


Since c(t) and d(t) are statistically independent, IP^(x) 
can be written 


oo x-X 

IP (x) = ; P, (X) / P (Y)dydX 

n \=— 00 y=i—co 


(A-26) 


Differentiation of equation (A-26) with respect to x gives 

Pu(x) 

00 

P (x) = / P^(X)P^(x-X)dX ( 

^ —00 

Substitution for from equation (A-20) gives 


P^(x) 


= / / Pa(X)[P^(Y)Pj^ 

X=-oo Y=0 


(Xzi) 


+ P^(-Y)Pj, 


(^)l 


^ dX 
Y 


(A-28) 


Now substitute the expressions for 

equations (A— 17) , (A— 18) , and (A— 19) . 

to X and make the substitutions 


, P^ , and P from 
' b c 

Integrate with respect 
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Y 


(A-29) 
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a /2 
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^ cr, 

a b 
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1/2 


+ ''d 


(A-30) 


(A-31) 


Where is the standard deviation of u(t) and R is the 

ratio of to 0 ^ . The final result is 


0 P 
u u 
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U+Rl) 


1/2 


7T 
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0 1+2 CR 
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exp[-C^-i (;^) \ ] 
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(A-32) 


2 '0 ^ 2 2 J 

u (1+2C R^) 


This result is presented in normalized form so that the choice 
rio influence on the resulting curve. Note the 

presence of R as a parameter of the equation and the absence 
of any dependence upon the filter transfer functions other than 
through R . Equation (A-29) is easily evaluated numerically, 
and results for various values of the standard deviation ratio 
R are presented in figure A-9, 

Note that R is the ratio of 0 ^ to 0 ^ and recall that 

d(t) is a Gaussian process while c(t) is the product of 
Gaussian processes. It is shown in reference A-5 that the 
probability density of c(t) is characterized by the modified 
Bessel function of second type and order zero, it is therefore 
called a Bessel function process. R is the ratio of the 
amounts” of these two processes used to generate u(t) . When 
R is very large, u(t) is composed almost entirely of the 
Bessel function process. Similarly, when R is very small 
u(t) IS composed almost entirely of the Gaussian process. 

The curves of figure A-9 reflect this result. 

2.) Cumulative Probability Distribution, 3P 

u 

The cumulative probability distribution of the system 
output can be obtained from equation (A-32) by integration. 
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(A-33) 


IP 


u 



x/a 


u 


•5 + ^0 


P 

u 


(^) dY 


Substitute for from equation (A-32) , interchange the order 

of integration, and make the substitution for Y 


X = Y 




(A-34) 


The result is 


1/2 


X r 1+R . 

CO a ^ _.2„2^ 


IP (JL) = .5+-/ 1+2C R exp(-X^)exp(-C )dXd? (A-35) 

u ^5=0 X=0 


Integrate with respect to 


and obtain the final result. 

2 1/2 


,1+R . 

CO X {— 2~) 

I> (-^) = .5 + — / exp (-5) erf[ o 2 1/2 


] d? (A-36) 


u 


u 


where erf denotes the error function. Equation (A-36) is 
easily evaluated by standard numerical techniques. Results 
for various values of R are given in figure A— 10. 

3.) Normalized Central Moments, . 

The normalized central moments of the system output can 
be obtained from equation (A-32) by integration. Since the 
average value of u(t) is zero, its nth normalized central 
moment is given by 


M 
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00 


f 


n 
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if-) dx 


(A-37) 
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Substitute for from equation (A-32) . Interchange the order 

of integration and integrate by parts to obtain 


M = 


n 

(n even) 


2[ (1) (3) •• • (n-1) 1 

/iT (1+R^)"/2 


00 n/2 

/ [2?^R^+1] exp(-5^) dC 


(A-38) 


M = 0 
n 

(n odd) 

Integrate equation (A-38) to obtain 


(A-39) 


M = 


n 

(n even) 


[ (1) (3) • • • (n-1) 1 


n/2 

S 

j=0 


n/2 n-2j 
( j ) R 


(1) (3) **• (n-2j-l)] 

(A-40) 


n/2 

where ( j ) denotes the binomial coefficient. The dependence 
of upon R IS illustrated for and Mg in figure A-11. 


^ guinroary of Result s for Individual Components .- Expressions 
the power spectral density and probability distribution 
of the suggested turbulence simulation components have been 
derxved and are presented as equations (A-18 ) , (A-19) , (A-32) 

(A-40) . The last three of these are also presented 
graphically in figures A-9, A-10, and A-11. 


It IS important to note that the shape of each component's 
power spectral density is determined only by the transfer func- 
tions of the filters used to generate that component. Similarly, 
the probability distribution is determined only by the standard 
deviation ratio R . 


it is not necessary to introduce cross correlations 
between the various components, the results presented above are 
sufficient to permit construction of a complete turbulence 
simulation. Each of the five components listed at the begin- 
ning of this appendix would be generated by a system of the 
type shown in figure A-8. The power spectrum and probability 
distribution of each would be determined by the above mentioned 
equations . 
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Correlation of Components 

The problem of correlating the components of the proposed 
turbulence simulation will now be considered. Two types of 
correlation are required, cross spectral densities and patch 
correlation. 

Cross Spectral Densities .- This section discusses a tech- 
nique by which a non-zero cross spectral density relating the 
u and w components of the turbulence simulation can be intro- 
duced. The method is taken from reference A-5 and consists 
essentially of generating these two time histories from statis- 
tically dependent white noise sources. It will be shown that 
this can be accomplished without altering the probability dis- 
tribution or power spectral density of either component. 

Consider the system of figure A-12. Note that the right 
hand side, enclosed by dashed lines, consists of two systems 
of the type shown in figure A- 8 . However, the Gaussian white 
noise sources have been replaced by the linear filters 

and noise sources r) . appearing on the left of figure A-12. 

It is the purpose of these and to produce statistically 

dependent white noise time histories. If a set of transfer 
functions can be found such that the time histories denoted by 
Nj^(t) are Gaussian white noise, then all statements made in 

the previous section on Statistical Analysis regarding the two 
systems on the right hand side of figure A-12 remain true. In 
particular, the probability distributions and power spectral 
densities are unchanged. However, because the noise sources 
Til , TI 3 f and rig are used in the generation of both u(t) 

and w(t) , these two time histories will now be correlated. It 
will first be shown that non-trivial choices of the filters 
Hj^ are possible, then an expression for the resultant cross 

spectral density will be derived. 

1.) Construction of Gaussian White Noise from Non-White Noise 
Processes 

Consider the block diagram of figure A-13 which is taken 
from the top left corner of figure A-12. In that figure, 
and TI 2 independent Gaussian white noise sources. Hj^ 

and H 2 are linear filters. Certain conditions must be placed 

on H^ and H 2 such that N^(t) will be Gaussian white noise. 
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Recall that a Gaussian process remains Gaussian when passed 
through a linear filter. Furthermore^ a sum of independent 
Gaussian processes is Gaussian. Hence N^(t) must be Gaussian. 

Next consider the power spectral density of Nj^(t) . Recall 

^e expression relating the power spectral densities of the 
input and output of a linear system and note that the noise 
T]^ and H 2 are independent. Then the expression for $ is 

^1 

= |H3^(i2 7rf)| ^Kj^+|H2(i2Trf) (A-41) 

where and K 2 are the power spectral constants of ri^ and 

^2 respectively. Now require Nj^(t) to be white noise so 

that 4> is a constant 
N, 


(A-42) 

For simplicity require to equal K 2 . 

= ^2 (A-43) 

Then (A-41) becomes 

\e^ (i2irf ) I 2+1 (i 2 TTf ) I 2 = ^ 

^rme history Nj^ (t) will be Gaussian white noise if equa— 

tions (A— 43) and (A— 44) are satisfied. Requiring each of the 
white noise sources of figure A-12 and each of the filter pairs 
Hi - H 2 through satisfy equations (A-43) and 

(A-44) will cause all of the N^^(t) to be Gaussian white noise. 
2.) Cross Spectral Density 


The cross spectral density of w(t) and u(t) will now be 
computed. This will be done by first finding an expression 

correlation of w and u , then Fourier transforming 
to obtain the cross spectral density. ^ 
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In the time domain the response time history of a linear 
system is given by the convolution of its impulse response func- 
tion with the input time history. In this report the convolution 
operator is denoted by an asterisk. Therefore the response of 
linear filter to the input is given by 


00 

{rii*hi)t= (t-X)hj^ (X)dX 


(A-45) 


Note that convolution is associative. 

From figure A-12, expressions for u(t) and w(t) are: 
w(t) = 

t t 

+ [(Tl3*h3*g2) + <n4*h4*g2) 1 [ (ri5*h5*g3) ^ 

^ ^ (A-47) 


u(t) = (Tii*h/g 4 ) + (n8*hg*g4) 

t ^ 

+ [(n3*hg*g5)^+ (nio*^10*^5^^^ ^2*^12%^^^ 

(A-48) 

Now consider the cross correlation of these two time histories 


C (t) = E{u(t)w(t+T) } 

uw 


(A-49) 


Sxibstitition of equations (A-47) and (A-48) into (A-49) 
leads to a lengthy expression containing many convolutions. 

This equation will not be reproduced here, but , 

be followed in reducing it to the final result will be expl . 

write each convolution involving a white noise 

same form as equation (A-45) . Then each n has a t appearing 
in its argiunent. Now interchange the order of integration and 
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apply the expected value operator directly to the integrand. 
The expected value operates only on the time variable which 
appears only in the arguments of the white noise functions 
Recall that 


r 0 if i j 
(t+A) rij (t+B) } = | 

Lk^ 6(A-B) if i = j (A-50) 

where A and B are arbitrary constants, <5 is the Dirac delta 
function, and is the power spectral constant of the white 

noise source . Finally, take the Fourier transform with 

respect to t , recalling the Fourier identities 


00 

X (t-A)y (X)dX} = X(f)Y(f) 

oo 

X (t+A) y (A)dA } = X(f)Y(f) 


(A-51) 

(A-52) 


where _ denotes complex conjugate. The final result for 
^uw^^^ " K^[H^(i2TTf)G^(i2Trf)H^(i2Trf)G^(i2TTf)] 


+ [H^ (i27Tf)G2 (^^Trf )Hg (i2Trf)G^ (i2Trf ) ] * 


[Hg (i2TTf )G3 (i27rf)Hj^^(i2TTf)Gg (i2TTf) ] (A-53) 


This is the required expression for . Note that of 

” filters present in the generating system of figure 
A- 12 only the odd indexed appear in equation (A-53) . These 

can be chosen subject only to the condition (A-44) 

Hence considerable generality has been obtained in the specifi- 
cation of a cross spectral density. It is of course quite 
difficult to find the appropriate H. such that $ matches 

some particular function which is of interest. This report 
follows the method of reference A-5 in selecting involving 
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arbitrary constants. The resulting expression for is 

not as general as equation (A-53) , but does provide a relatively 
simple cross spectral shape with continuously adjustable param- 
eters. This approach to the problem of cross spectra was com- 
pletely adequate for the flight simulator experiment described 
in this report. 

Patch Correlation.- In addition to the conventional cross 
spectral density discussed in the preceding section, another 
type of correlation called "patch correlation" must be considered. 
If the simulated turbulence is to have a patchy character , 
observed time histories show that when a patch is encountered 
the increased activity should appear in all of the turbulence 
components at the same time. In an effort to obtain this sort 
of correlation without altering the probability distribution, 
power spectral density, or cross spectral density of the com- 
ponents, the following approach was used. Only the u and w 
gust components will be discussed here, but a similar method 
can be applied to all components of the simulated turbulence. 
Unfortunately, this technique is only useful in the case of a 
digitally generated simulation. 

Figure A-14 shows the cross correlated u and w gust com- 
ponents of the turbulence generator as in figure A-12, but 
with the noise sources Hg / > and t 1 j ^2 replaced by signal 

paths from 1 I 2 / t and rig • As before, the are inde- 

pendent Gaussian white noise sources. In a digital system 
these noise sources produce not continuous functions of time 
but uncorrelated, normally distributed random numbers. The 
± 1 operators represent random sign functions which take on 
the gain values plus unity or minus unity as each sample of 
turbulence is computed. Since successive samples of a white 
noise process are uncorrelated, these random sign functions 
alter neither the Gaussian nor the white noise nature of the 
inputs to Hg , H^q , and ‘ Furthermore, if and 

denote respectively the input and output of a random sign 
function, then r\^ and are independent in the sense that 

E{n^Tly} = E{n^}E{Tiy} = 0 (A-54) 

However, they are correlated in the sense that 

= E{Tiy^} (A-55) 
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implies that the cross spectral density of 
u(t) and w(t) remains unchanged by this patch correlation 
technique. However, in view of equation (A-55) , it is possible 
toat large excursions of w(t) will tend to be acconpaSirby 
excursions of u(t) , This has not been proven, but 
examination of sample time histories produced by the proposed 
turbulence simulation reveals that patches of activity do tend 
to appear on all components at the same time. it will be shSwn 
later how the above described procedure was extended to include 
all components of the simulation. 


Rotation of Coordinates to Provide Cross Spectra Relating 
Components Other Than u and w 

+■>, flight simulator experiment described in 

seve^rar??fi°^ report was to compare the realism of 

i turbulence models with that of actual atmos- 

f artificial models were to differ from 
their probability distributions, 
e efore it was necessary to eliminate as many other differ- 

possible. One of these was the cross spectra of the 
turbulence components. 

cross spectra has already been discussed in 

made^irthat However, the remarks 

A referred to the gust components measured in 

system. Unfortunately, the real tur- 

were^all used in the flight simulator experiment 

TS obtain^ in a vehicle fixed coordinate system. Since 

^elaJivS any heading 

relative to the mean wind direction, non-zero cross spectra 

a?Z^po2Iible^°''®"^® longitudinal and vertical 

There are two ways to match the cross spectra of the arti- 
ficial models to those of the real time histories: 

simulated vehicle at various headings relative 
essentially modeling the manner in which 
the real turbulence data were obtained; 

(2) rotate the coordinates in which the artificial models 
coL?an?!^^®'^ simulated vehicle heading remains 

methods would simplify the turbulence genera- 
iZ v! u- '^ould require different vehicle headings^for 

Doini-^Z?® tested. This is undesirable from the^stand- 

point of providing a uniform flight task for all models. 


no 



The second approach complicates the turbulence simulation 
but permits the seune flight task to be used for evaluating 
each model. It is therefore the method adopted in this report. 

Equations describing the influence of a coordinate rotation 
in the u - v plane upon the spectra and cross spectra observed 
by the simulated vehicle will now be derived. Figure A-15 
shows the geometry of the rotation. The gust components in 
the vehicle coordinate system are; 


u^ = u cos 

+ V 

sin 

'p 

(A-56) 

v' = V cos 

- u 

sin 


(A-57) 

w 

= w 



(A-58) 


The spectra and cross spectra of these components are: 

2 . 2 

$ $ cos + $ sin iIj 

u u V 

2 2 

$ ^ $ cos sin 

V V u 


$,=(!) 
w w 


t = tt ) sini241 

u V V u 2 


(A-59) 

(A-60) 

(A-61) 


(A-62) 


$ ^ ^ $ cos (A-63) 

u w uw 


$ ^ ^ =-4 sin \l> (A-64) 

V w uw ^ 

Using equations (A-59) through (A-64) it is possible to work 
backward from the known spectra and cross spectra of the real 
turbulence time histories in vehicle coordinates to a suitable 
set of spectra and a u - w cross spectra in earth-fixed coor- 
dinates. Exact solutions are generally impossible because the 
six equations involve only five unknowns/ but reasonable results 
may be obtained for certain spectra shapes, and this was possible 
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for the time histories used in this report. It should be empha- 
sized that this procedure was necessary only because the objec- 
tive of the flight simulation was to evaluate the realism of 
various turbulence models. The coordinate rotation used here 
would not be necessary in a typical flight simulator experiment. 


System Diagram and Specialization to Produce 
Specific Turbulence Models 

The equations and techniques of the preceding sections 
have been applied to produce the complete system of figure A-16 
and Table A-2. The results of reference A-5 have been used 
extensively. 

The system shown is suitable only for digital generation 
techniques because of the random sign functions required by 
patch correlation. If these are eliminated the system could 
be used with an analog computer, but additional noise generators 
would be required in order to prevent unwanted cross correlations. 

The coordinate rotation applied to the u and v gust 
components follows the method suggested in the previous section 
on "Rotation of Coordinates". This rotation is necessary only 
for the specific experiment described in this report and is 
not required for a typical flight simulation. 

The proposed turbulence simulation can be simplified greatly 
if the u-w cross correlation is not required. In this case 
only the filters enclosed by the dashed lines in figure A-16 
are required. 

The five output time histories of the complete system have 
the frequency and probability characteristics described pre- 
viously in this appendix. However the scale lengths, cross 
correlation parameters, and coordinate rotation angles required 
to model the real time histories used in the flight simulator 
experiment have not been specified. These pareuneters are 
tabulated in Table A- 3 for the two LO-LOCAT time histories 
used in this report. 
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TABLE A-2.- TRANSFER FUNCTIONS OF FIGURE A-16 
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'=3'V^6> 
2L s 

(1 + — ^) 


(0.004842) (1+0. 7425s) (1+0. 2270s) 
4^(1+R^) (1+1. 2569 s) (1+0. 1341s) (1+0.2653s)^ 


G23(s) = 


(.016219) 


(K 3 +K^) d+R"^) (1+4. 3993s) 


(1+0. 2391s) (1+0. 9324s) 

(Ks+Kg) (1+1. 7596s) (1+0. 1383s) (1+0. 2823s) 


o (0.006331) (1+0. 7314s) (1+0. 1154s) 


'Kj^( 1+R'^) (1+1. 2410s) (1+0. 06803s) (1+. 19733s) 


^ 26 ^") = 


a^R/C^ (.021162) 


(K^+K^) (1+R^) (1+4. 3993s) 


(1+0. 1171s) (1+0. 9166s) 

:^(K^+K^) (1+1. 7286s) (1+0. 0691s) (1+0. 2065s) 
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TABLE A-3.- FILTER COEFFICIENTS FOR SPECIFIC TURBULENCE 
SIMULATIONS 


A. are arbitrary scaling factors, may be set equal to unity 
^ if desired 

K. must be set equal to power spectral constant of correspond- 
^ ing white noise source 

u^ = 36.02 m/sec (118 ft/sec) in all cases 

= .765 m/sec (2.51 ft/sec) in all cases 

0 ^ = .832 m/sec (2.73 ft/sec) in all cases 

0 ^ = .579 m/sec (1.90 ft/sec) in all cases 


Models Based on Real Turbulence Scimple from LO-LOCAT Test 1140, 
Condition 21, Category 412112 



GAUSSIAN 

MATCHED 

NON-GAUSSIAN 

MORE PATCHY 
NON-GAUSSIAN 


170.7 m (560 ft) 

170.7 m (560 ft) 

170.7 m (560 ft) 

\ 

141.7 m (465 ft) 

141.7 m (465 ft) 

141.7 m (465 ft) 

K 

141.7 m (465 ft) 

141.7 m (465 ft) 

141.7 m (465 ft) 


1.9 

1.7 

1.6 

c„ 

_ _ _ 

3.6 

3.1 

K 




ip (deg.) 

45 

45 

45 


Models Based on Real Turbulence Sample from LO-LOCAT Test 1141, 
Condition 21, Category 413212 



GAUSSIAN 

MATCHED 

NON-GAUSSIAN 

MORE PATCHY 
NON-GAUSSIAN 

^u 

^w 

^G 

'Jj (deg.) 

90.2 m (296 ft) 
108.2 m (355 ft) 
36.0 m (118 ft) 
2.9 

1 

0 

90.2 m (296 ft) 
108.2 m (355 ft) 
36.0 m (118 ft) 
2.55 
5.05 
0 

90.2 m (296 ft) 
108.2 m (355 ft) 
36.0 m (118 ft) 
2.55 
5.05 

0 
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/ DATA FROM REF. A-4 
'1 VERTICAL GUSTS 

KANSAS PLAINS 
STABLE CONDITIONS 

^76 M. ALTITUDE 

2Z9 M. ALTITUDE 

1 1 1 1 1 1 1 

lA 

/ 
y J 


-6 -4 -Z 0 2 4 6 


NORMALIZED GUST VELOCITY.;^ 

cr^ 


FIGURE fl-1. TYPICAL CUMULATIVE PROBABILITY 
DISTRIBUTIONS OF ATMOSPHERIC TURBULENCE 


119 




FIGURE fl-2. TYPICAL PROBABILITY DENSITY 
DISTRIBUTIONS OF ATMOSPHERIC TURBULENCE 
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FIGURE fl-3. PATCHY CHARACTERISTICS OF HDOIFIED BESSEL 
FUNCTION MODEL COMPARED WITH THOSE PRODUCED BY THE 
GAUSSIAN MODEL AND ACTUAL ATMOSPHERIC TURBULENCE 



GAUSSIAN 

WHITE NOISE LINEAR 
SOURCES FILTER 


OUTPUT 

u(t) 



FIGURE A-4. LINEARLY FILTERED 
GAUSSIAN WHITE NOISE TURBULENCE 
SIMULATION 



OUTPUT 

u(t) 


FIGURE A -5. A NON- GAUSSIAN SIMULATION 
BASED UPON THE MULTIPLICATION OF TWO 
GAUSSIAN PROCESSES 
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CUMULfiTIVE PROBABILITY. P( 



FIGURE A-6. COMPARISON OF GAUSSIAN AND 
MODIFIED BESSEL FUNCTION CUMULATIVE 
PROBABILITY DISTRIBUTIONS 
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FIGURE fl-7. COMPARISON OF GAUSSIAN 
AND MODIFIED BESSEL FUNCTION 
PROBABILITY DENSITIES 
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GAUSSIAN 
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FIGURE A-8. PROPOSED NON-GAUSSIAN TURBULENCE SIMULATION 




3 


L 2 

NORMALIZED GUST VELOCITY, 


W 

CTw 


FIGURE fl-9. COMPARISON OF PROBflGILITY 
DENSITIES FOR VARIOUS VALUES OF THE 
STANDARD DEVIATION RATIO 




FIGURE fl-10. COMPARISON OF CUMULATIVE 
PROBABILITY DISTRIBUTIONS FOR VARIOUS 
STANDARD DEVIATION RATIOS 
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MOMENT VALUE 



STflNDPRD DEVIATION RRTIO. R 


FIGURE fl-11. DEPENDENCE OF THE FOURTH 
RND SIXTH NORMALIZED CENTRAL MOMENTS 
UPON THE STANDARD DEVIATION RATIO 
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FIGURE A- 13. SYSTEM FOR 
THE GENERATION OF 
GAUSSIAN WHITE NOISE 
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FIGURE A- 14. SYSTEM FOR THE 
GENERATION OF CORRELATED PATCHES 
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MEAN 



FIGURE A- 15. COORDINATE ROTATION 

GEOMETRY 
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FIGURE A- 16 BLOCK DIAGRAM OF COMPLETE 

SIMULATION 
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APPENDIX B 


EQUATIONS OF MOTION 


The following equations of motion are essentially those 
provided by NASA Ames for representing the Twin Otter for use 
on their simulators. They have been modified slightly to 
include the turbulence inputs. 

The total (non-perturbation) forces and moments on the air- 
craft due to aerodynamic and engine forces are: 


F 

Tx 

= qSC^ t Fg 

X 

(B-1) 

Ty 

= 5sc^ + 

(B-2) 

F 

Tz 


(B-3) 

Tx 

= isbC^ + 

(B-4) 

^Ty 

= qScC + M 

m Ey 

(B-5) 

^z 

= qSbC +M„ 

^ n Ez 

(B-6) 


The axes system (x,y,z) is a stability-axes system with the 
usual convention for signs. The dynamic pressure is computed 
using 

where (Ug , Vg , Wg) are the body-axes vehicle velocities with 

respect to the gust air. Thus q contains the effects of gusts 
as well as changes in airspeed. 

The engine forces and moments are 




(B-7) 
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(B- 8 ) 


F = -qSC C (6T-6T) 
EZ L^ip 


^Ev = (<5T-6T) 

"'CT ^6t 


(B-9) 


=0 

Ey Ex Ez 


(B-10) 


The thrust sensitivity is 


= 6.5617 + 0.0021) (B-11) 


qs " V 




The aerodynamic force and moment coefficients are given by 


Cx = - ^ 


C = C 6 + 


'y yg 2v^w yr ^ 


(Cy r.Cy P) .Cy^^aa 


+ C «r + ^ C„ (v„ - V ) 
^Sr ''rW ^te 9c 


c = - C, +C, (C o+C q) 


+ C 6 e + C (w^ - w ) 
^6e '^RW ^ta ^t ^c 


+ — ^ C, . (w„ - w ) 

2vZ 7 ^ta ^t ^c 
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3 ^''rW ""r ^6a 


(C^ r + C^ p) + c 6a 


6r ug wg 


(B-16) 


C_ = C_ a + 

V ^^RW “‘a “‘q 


m -m. - 2V:;:: «e 

6e 


+ 


(u - u ) + c (w - w i ^ 

m. a. a ' m f / 


'^RW ■ "'tu' ^t ^c' "^c 


+ 0 . (w - w ) 

2 VrW ta 


(B-17) 


^ ^ 2C:: P> C_ 6a 


n n 


3 


RW “r 


n 


6a 


^ ^n '*' V^ ^r, ~ ) + C (t) 

«r ''rw "t6 


+ C^ (t) 
wg 


(B-18) 


Here, a and 3 
gusts; 


dsfined so as to incXude the effects of 


w 

a = arc tan — 
^B 


(B-19) 


3 = arc tan 


B 




(B-20) 


Th6 forces and moments (B-1) to (B-6) were fed into a NASA 
digital program which broke up the forces into a north- east -down 
axes system, added in the effects of gravity and earth's rota- 
tion, and performed the necessary integrations to solve the 
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equations of motion. The equations used in this digital program 
were all exact as all nonlinearities were included. The details 
of these equations are relatively straightforward and will not 
be gone into here. 
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APPENDIX C 


MOTION COMPENSATION 


Three different types of motion compensation were used; 
acceleration limits, a "quickening" circuit to improve the hiqh- 
frequency response of the simulator, and a washout circuit to 
filter out low-frequency motions. The arrangement of these cir- 
cuits IS shown in figure C-1. 


The drive-system for the simulator was a position 
i.e., a position was commanded by the computer and the 
drive-system tried to match the commanded position. 


-drive, 

simulator 


The inputs to the motion-compensation network were the 
computed Euler-angle accelerations V / V , and V , and the 
computer accelerations , and Vj^ in the north-east- 

down axis system. The accelerations / V^. , and V were 

computed by taking the total forces (in body-axes) on the air- 

them into forces in the north-east-down 
axis system by using the washed-out Euler angles of the simula- 
'J’w ' ®w ' K ' dividing by the airplane mass m . 


Acceleration Limiting 

In order to prevent the drive system from shutting down 
because the commanded accelerations were too large, limits were 
placed on the acceleration signals. Because the simulator was 
ariven by position rather than by acceleration, the limiting 

was done on the computed values of , V^. , and Vj^ instead 

of on actual- cab accelerations in order to reduce the complexity 
Thus, because the signals are passed through a quickening cir- 
cuit before reaching the simulator drive, the commanded accel- 
erations sometimes exceeded the acceleration limits set. The 
Imiting was done digitally (as was all the compensation) , and 
the following values were used for limits; 


= ±1.524 m/sec^ 
Vg = ±1.524 m/sec^ 
Vjj = ±0.9144 m/sec^ 
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Quickening 


A quickening, or lead, circuit was used on all six signals 
in order to improve the high-frequency response of the simulator. 
The accelerations were first integrated twice to give distance 
and angle commands. The quickening input/output relations were 
of the form 


'out 


= e. + 
in 


K, 


e. + 
in 


K 


DD 


e . 
in 


Taking the Laplace transform. 


^ = 1 + K s + Kdd 

in 


Note that this gives increasingly higher gain as frequency is 
increased, and is partly the cause of the excessive high-fre- 
quency motion that the pilots complained about in the simulation. 
The following values were used for and in each of 

the channels. 


Channel 


^DD 

4> 

0.08 

0. 

0 

0.18 

0.012 

t 

0. 

0. 

*NORTH (longitudinal) 

0.16 

0.04 

^EAST (side-motions) 

0.14 

0.037 

*DOWN (''ertical) 

0.25 

0. 


Washout 

Washout was necessary in order to prevent a low-frequency 
acceleration from building up motions that would cause the simu- 
lator to hit its travel limits. The washout used had an input/ 
output Laplace transform of the form 


^in s^+As + B 
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This provided unity gain at the high frequencies, and a gain 

goxng to zero as frequency went to zero. The values of A and B 
used were as follows. 
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OISTflNCES DISTANCES INPUTS TO 

M] ANGLES AND ANGLES SIMULATOR 
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FIGURE C-1. DETRILS OF MOTION COMPENSATION USED IN SIMULATION 






APPENDIX D 


PILOT QUESTION SHEET 


Twin Otter 6-Degree Simulation, June-July 1971 


Flight Number; 
Pilot; 


Date ; 


1 . Turbulence Intensity ; 


Light 
( ) 


Moderate 

( ) 


2. Realism of Turbulence; 


Very Good 

{ ) 


Good 

( ) 


Severe 
( ) 


Fair 

( ) 


Extreme 

( ) 


Poor 

( ) 


Very Poor 

( ) 


3. Correctness of Relative Amplitudes of Disturbances ; 


Roll 

Pitch 

Yaw 

Heave 

Side-Force 


Not Enough About Right Too Much No Comment 

; ! ( ) ( ) ( ) 

( ) ( ) ( ) 

( ) ( ) ( ) 

( ) ( ) ( ) 

( ) ( ) ( ) ( ) 


4 • Characteristics (Variations of Intensity— Burst ) 


Much too 
Continuous 
(monotonous) 
( ) 


A Little too About A Little No 

Continuous Right too Patchy Comment 

( ) ( ) ( ) ( ) 


5. Frequency Content of Turbulence ; 

Not Enough About Right Too Much No Comment 
Low Frequencies ( ) ( ) ( ) ( j 

High Frequencies ( ) ( ) ( ) ( j 

6. Other Comments about Realism (of turbulence or aircraft 
simulation) ; 


7. Pilot Estimate of Work Load: 


Very Easy 

( ) 


Easy 

( ) 


Average 

( ) 


Difficult 

( ) 


Very Difficult 

( ) 
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8. Pilot Estimate of Task Performance (integral-squared error 
for ILS tracking task) : 

Very Good Good Average Poor Very Poor 

( ) ( ) ( ) ( ) ( ) 

9. Cooper Rating (rate aircraft as flown under these conditions 
of turbulence) : 
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APPENDIX E 


OTHER RESULTS OF PILOT QUESTIONNAIRE 


The following plots show the pilot answers to the question 
sheet which were not discussed in the main text. They do not 
show any consistent differences between the various turbulence 
models and are shown here only for completeness. Also shown 
are plots of the ras values of the angles <(>, 6, and . Again, 
there are no consistent differences between the rms values from 
one model to the next. 
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PILOT RRTING -o PILOT RRTING 


PILOT DflTfl POINTS 
X PILOT fl 
□ PILOT B 


PILOT fl PILOT B 




+95 X CONFIDENCE LIMIT FOR MEAN 

MEAN OF DflTfl POINTS 

-95 X CONFIDENCE LIMIT FOR MEAN 


NO COMMENT 

TOO MUCH 

ABOUT 

RIGHT 

NOT ENOUGH 


NON-GflUSSIflN NON-GflUSSIflN 

MODEL 



GAUSSIAN REAL MATCHED MORE PATCHY 


FIGURE E-1. CORRECTNESS OF PITCH 


ILOT DATA POINTS 


X PILOT A 


0 PILOT B 


PILOT fl PILOT B 




+95 X CONFIDENCE LIMIT FOR MEAN 

MEAN OF DflTfl POINTS 

-95 X CONFIDENCE LIMIT FOR MEAN 


NO COMMENT 

TOO MUCH 

ABOUT 

RIGHT 

NOT ENOUGH 


NON-GflUSSIflN NON-GflUSSIflN 

MODEL 



GAUSSIAN REAL MATCHED MORE PATCHY 


FIGURE E-2. CORRECTNESS OF HERVE 
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PILOT RATING -o PILOT RATING 


PILOT OflTfl POINTS 
X PILOT fl 
□ PILOT B 


PILOT fl PILOT B 




+95 X CONFIDENCE UNIT FOR HERN 

NEON OF OOTO POINTS 

-95 X CONFIDENCE UNIT FOR NEON 


NO CONNENT 

TOO NUCH 

ABOUT 

RIGHT 

NOT ENOUGH 


NON-GOUSSION NON-GOUSSION 

MODEL 



N/V 



vrv/ 






V/ 


kk 




> 

< 

V/ 

> 

''m 

J- 



/ 



— 

alt- 


— 

□ j- 

Timrn 1 








GAUSSIAN REAL NOTCHED NORE PATCHY 


■IGURE E-3. 


CORRECTNESS OF SIDE-FORCE 


ILOT DATA POINTS 
X PILOT A 
□ PILOT B 


PILOT A 




PILOT B 

+95 X CONFIDENCE UNIT FOR NEON 
NEON OF DATA POINTS 
95 X CONFIDENCE LINIT FOR NEON 


VERY DIFF 
DIFFICULT 
AVERAGE 
EASY 
VERY EASY 







■ 1 X 
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